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1,  INTRODUCTION 


The  concent  of  ae • r.g  the  electric  field  of  Ihp  c.  ..o- 
sphere  for  terrain  avoidance  had  its  origin  in  me  sucrcstful 
use  of  the  atmospheric,  ph.'tric  field  to  stabilize  the  flight 
of  model  airplanes  by  M.  L.  Hill  (Ref.  1).  Using  r.labvely 
simple  electronics  and  servo-  * ichanisr.v.,  Hi)/  snowed 
that  it  was  possible  to  fly  a modt1  a:rc.af>  ...  remotely  pi- 
loted vehicle  (RPV)  along  *he  eauipotential  surfaces  that 
exist  in  the  atmosphere  parallel  to  the  earth's  surface. 

This  "electrostatic  stabilization"  of  flight  over  relatively 
level  terrain  has  been  demonstrated  to  be  possible  wher- 
ever the  so-called  "fair  weather"  electric  atmospheric  field 
exists  — which  is  about  90%  of  the  time.  Since  the  conduc- 
tivity of  the  ground  (or  any  orographic  protrusion)  is  so 
much  greater  than  that  of  the  atmosphere,  the  boundary 
condition  for  Laplace's,  or  more  properly  Poisson's  equa- 
tion describing  the  atmospheric  potential,  is  that  the  poten- 
tial at  the  earth's  surface  and  the  surface  of  any  protrusion 
is  a constant  conventionally  taken  as  zero.  Thus  at  the 
earth's  surface,  including  all  protrusions,  whether  man- 
made or  natural,  there  exists  an  equipotential  surface,  and 
the  atmospheric  field  in  the  neighborhood  of  these  protru- 
sions will  be  distorted.  The  postulation  of  the  use  of  the 
atmospheric  field  for  terrain  avoidance  depends  upon  the 
detection  and  use  of  these  spatial  disturbances. 


The  first  problem  in  this  investigation  was  to  de- 
termine the  magnitude  and  direction  of  the  changes  in  the 
atmospheric  electric  field  in  the  space  surrounding  the  oro- 
graphic protrusion.  Later  in  this  report,  it  will  be  shown 
that  it  is  possible  to  obtain  computer  solutions  and  in  some 
cases  analytical  solutions  of  the  atmospheric  potential  field 
in  the  neighborhood  of  realistic  orographic  protrusions. 

Ref.  1.  M.  L.  Hill,  "Introducing  the  Electrostatic 
Autopilot.  " Astronautics  and  Aeronautics.  22-31  November 
1S72. 
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provided  certain  assumptions  are  made.  These  compu- 
tations make  it  possible  to  specify  the  changes  relative  to 
the  unperturbed  atmospheric  field  in  the  vertical  and  hori- 
zontal potential  gradients  as  a function  of  the  horizontal 
distance  from  the  protrusion  at  a specified  altitude. 

Given  the  results  of  these  computations  two  questions 
exist  that  can  only  be  answered  by  experiment.  Are  the  as- 
sumptions made  for  the  computations  realistic  enough  to 
yield  results  applicable  in  the  real  world?  And  if  the  po- 
tential gradients  are  found  to  be  as  predicted  by  Poisson's 
equation  is  it  possible  to  devise  reliable  instrumentation 
that  will  measure  them  accurately?  For  the  limited  range 
of  parameters  measured,  this  report  will  show  that  the 
answers  are  affirmative.  The  limitations  in  the  experi- 
mental work  reported  are  not  inherent  in  the  problem,  but 
reflect  a limited  effort. 


1 


The  emphasis  of  the  present  investigation  has  been 
on  the  atmospheric  electric  field  at  a relatively  low  altitude, 
i.  e. , below  150  meters.  Measurements  below  3 meters 
were  made  in  the  vicinity  of  the  Applied  Physics  Labora- 
tory in  Howard  County,  Maryland.  Measurements  at  higher 
altitudes  were  made  of  the  South  Mountain  area  in  western 
Maryland  from  a Cessna  337  aircraft. 

In  general,  the  experimental  program  has  been  lim- 
ited to  the  "fair  weather"  atmospheric  electric  field.  In 
the  study  of  atmospheric  electricity  the  term  "fair  weath- 
er" usually  refers  to  those  occasions  when  the  electric  field 
vector  is  directed  toward  the  earth  and  there  is  no  precipi- 
tation or  electrical  storm  activity  in  the  area.  There  is  a 
high  degree  of  correlation  between  the  fair  weather  electric 
field  and  "fair  weather"  in  the  meteorological  sense,  but 
the  atmospheric  electric  field  does  occasionally  change  sign 
during  meteorologically  fair  weather.  As  is  shown  later 
the  fair  weather  atmospheric  field  is  subject  to  significant 
fluctuations  both  in  time  and  space.  Although  field  reversals 
and  abnormally  high  field  strengths  are  usually  associated 
with  thunderstorm  activity  or  the  passage  of  storm  fronts, 
these  conditions  have  been  avoided  for  reasons  of  experi- 
mental convenience.  Even  under  the  most  favorable 
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conditions  the  atmospheric  electric  field  is  difficult  to  mea- 
sure as  the  atmospheric  potential  is  regarded  as  a very 
high  impedance  voltage  source.  Also,  reliable  measure- 
ments are  quite  difficult  to  make  during  periods  of  precipi- 
tation when  insulation  and  triboelectric  charging  problems 
become  severe.  Measurements  during  periods  of  lightning 
activity  require  instrumentation  of  enormous  dynamic  range. 
This  instrumentation  must  be  protected  against  being  de- 
stroyed by  the  high  electrical  field  and  current. 

For  the  purposes  of  this  investigation  this  restriction 
to  the  fair  weather  electric  field  is  believed  to  be  realistic 
since  the  application  of  the  atmospheric  electric  field  for 
stabilization  or  terrain  avoidance  will  be  much  more  diffi- 
cult or  even  impossible  under  other  atmospheric  conditions. 
The  percentage  of  time  during  which  the  fair  weather  field 
exists  and  the  geographical  extent  over  which  the  electric 
field  is  uniform  enough  to  be  useful  are  crucial  questions 
when  requirements  for  the  mission  of  an  aircraft  or  missile 
are  involved.  The  sensitivity  of  the  instrumentation  re- 
quired for  the  detection  of  the  atmospheric  electric  field 
perturbations  for  purposes  of  terrain  avoidance  is  a function 
of  the  dynamic  capabilities  of  the  aircraft  or  missile  and 
the  nature  of  the  protrusion.  These  capabilities  are  not  ex- 
plicitly considered  in  this  report. 
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Although  the  importance  of  electrical  phenomena  in 
the  atmosphere  was  appreciated  in  the  eighteenth  century, 
and  some  significant  investigations  were  made  in  the  nine- 
teenth century,  the  present  understanding  of  the  atmo- 
spheric electric  field  dates  from  about  1920  (Refs.  2 and  3). 
At  about  50  km  altitude,  the  atmosphere  is  (in  effect)  a per- 
fect conductor;  this  conducting  layer  (the  electrosphere)  to- 
gether with  the  earth  forms  a spherical  condenser.  The 
potential  of  the  electrosphere  is  about  300  000  volts  posi- 
tive relative  to  the  earth  which  is  conventionally  taken  to 
be  at  a potential  of  zero.  Therefore  an  electric  field  exists 
in  the  atmosphere  that  in  any  natural  sysi’m  of  coordinates 
would  be  negative  (i.  e. , pointed  toward  the  earth)  but  which 
is  usually  denoted  as  "positive”  in  the  study  of  atmospheric 
electricity.  This  confusion  of  sign  can  be  avoided  by  re- 
ferring to  the  atmospheric  electric  field  as  the  "potential 
gradient.  " So  designated,  the  atmospheric  potential  gradient 
is  for  the  normal  fair  weather  field  positive  in  conformity 
with  the  usual  conventions  cf  electrical  theory. 

The  flow  of  current  between  the  electrosphere  and 
the  surface  of  the  earth  is  carried  by  positive  ions  moving 
downward  and  negative  ions  moving  upward.  This  current 
density  is  found  to  be  constant  or  nearly  constant  as  a func- 
tion of  altitude  over  any  given  area  measured  at  a given 
time.  The  value  of  this  current  density  is  between  1 * 1 0~  1 2 
and  4 x 10-12  A/m2. 


Ref.  2.  J.  A.  Chalmers,  Atmospheric  Electricity. 
Pergamon  Press,  New  York,  K.  Y. , 1967. 

Ref.  3.  H.  Israel,  Atmospheric  Electricity.  Vols.  I 
and  II,  U.  S.  Department  of  Commerce,  National  Technical 
Information  Service,  Springfield,  Va.,  1973. 
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The  conductivity  ol  the  atmosphere  increases  with 
altitude.  This  conductivity  results  in  a decrease  in  the  ab- 
solute value  of  the  electric  field  with  altitude  since 


J=Xe,  (1) 

where  X is  the  conductivity.  J is  essentially  a constant. 
ience  E varies  inversely  with  X.  The  value  of  X near  the 
earth  is  typically  2 x 10"14  mho/m  so  that  for  J = 2 x 10-12 
A/m2  the  result  is 


E = ^ = 100  V/m  . (2) 

To  avoid  the  sign  confusion  inherent  in  atmospheric  electric- 
ity, where  altitude  is  taken  as  a positive  coordinate,  it  is 
customary  to  write  Eq.  (2)  as 


F = i, 


where  F is  the  potential  gradient  and  I is  a positive  current 
density  flowing  toward  the  earth. 

An  obvious  consequence  of  the  decrease  in  ihe  po- 
tential gradient  with  altitude  is  that  the  electric  field  has  a 

ip 

divergence,  (i.  e. , f 0),  and  from 


7-E  = o/Co,  (4) 

where  o is  the  charge  density,  it  is  clear  that  the  atmo- 
sphere is  not  electrically  neutral  but  must  under  normal 
conditions  contain  a net  positive  charge  per  unit  volume. 

It  is  clear  that  the  atmosphere  is  a conducting 
medium  of  very  high  resistivity.  For  aim2  column  of  the 
atmosphere  extending  from  the  earth's  surface  to  the  electro- 
sphere, the  columnar  resistance  is  of  the  order  of  101? 
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ohms.  Most  of  this  resistance  lies  in  the  lower  atmosphere, 
approximately  one  half  of  it  below  3000  meters.  A repre- 
sentative value  of  the  resistivity  of  the  atmosphere  would 
be  about  5 x 10*3  ohms/m  near  the  earth's  surface. 

There  is  a great  variability  in  the  atmospheric  con- 
ductivity, and  no  good  general  law  for  its  variation  with 
altitude  exists  (e.  g. , see  Chapter  7 of  Ref.  2).  A common 
method  of  fitting  empirical  data  has  been  to  use  a series  of 
exponentials,  of  which  the  first  term  predominates  at  lower 
altitudes.  A typical  equation,  to  the  first  power  in  altitude 
(z),  would  be 

X = 2.  5 x 10‘14  exp  (0.  00025z)  , (5) 

which  provides  a useful  fit  to  the  data  for  the  first  few  thou- 
sand meters  of  altitude. 

Since  in  Eq.  (3)  the  current^ density  is  the  relatively 
constant  term,  the  fluctuations  in  F are  the  consequences  of 
the  fluctuations  inX.  In  the  lower  atmosphere  electron  life- 
time is  less  than  a microsecond,  and  the  charge  carriers 
are  positive  and  negative  ions.  The  conductivity  is  given  by 

X = X+  + X_  = n+c^+  + n_«H_  , 

where  n+  and  n are  the  concentrations  of  positive  and  neg- 
ative ions  and  M+  and  I*,  are  the  mobilities.  More  properly 
since  there  is  sufficient  reason  to  suppose  that  several 
species  of  ions  exist  at  any  one  time,  the  conductivity  should 
be  written  as 

X = S e(n+M+)  +Ee(n_fO  , (6) 

i i 1 " j 

where  there  are  i species  of  positive  ions  and  j species 
of  negative  ions.  Thus  the  atmospheric  conductivity  in- 
volves two  questions:  first  how  is  ionization  produced  and 
second  what  are  the  properties  of  ions. 
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At  lower  altitudes,  the  contribution  to  ionization  by 
UV  radiation  is  negligible.  Over  the  ocean  surfaces,  the 
only  source  of  ionization  is  cosmic  rays.  Over  land,  ioni- 
zation is  produced  by  ground  radioactivity  and  by  radio- 
active gases  (principally  radon)  that  escape  from  the  ground. 
The  ionization  produced  by  ground  and  airborne  radioactivity 
is  very  variable,  but  is  typica’ly  several  times  that  pro- 
duced by  cosmic  rays  in  the  lower  atmosphere.  According 
to  Kef.  3 for  the  continental  atmosphere  the  rate  of  ioni- 
zation decreases  up  to  an  altitude  of  1 km  because  of  the 
decrease  in  earth-originated  radioactivity,  and  increases 
above  that  altitude  because  of  the  increased  effectiveness  of 
cosmic  radiation.  Since  the  rate  of  ionization  over  land  is 
5 to  10  times  that  over  the  ocean,  where  only  cosmic  radi- 
ation is  effective,  it  might  be  expected  that  the  atmospheric 
conductivity  would  be  much  greater  over  continental  areas. 

But  this  is  not  the  case  (in  fact  the  conductivities  over  land 
and  sea  are  approximately  equal)  because  of  the  complex 
nature  of  the  atmospheric  ions  and  the  interactions  between 
them.  Under  the  simple  fair  weather  assumptions,  where 
the  electrosphere  bears  a relatively  constant  potential  and 
the  conduction  current  to  the  earth  is  also  fairly  constant, 
it  is  the  wide  fluctuations  of  the  conductivity  (caused  by 
ionic  phenomena)  that  causes  the  fluctuations  in  atmospheric 
potential  (or  potential  gradient)  that  makes  terrain  avoidance 
or  any  other  application  of  the  earth's  atmospheric  field 
difficult. 

There  are  several  recognized  classes  of  ions  in  the 
atmosphere.  The  so-called  "small"  ions  are  known  to  be 
more  complex  than  a single  N2+  or  Oj*  molecule,  and  are 
thought  to  consist  of  from  four  to  ten  molecules  bound  in  a 
cluster  by  the  polarization  forces  exerted  by  the  initially 
single  charged  simple  molecule.  The  mobilities  of  these 
small  ions  are  typically  1 to  2 x lo-4  m2/V-s  or  1 to  2 
cm2/V-s  in  the  conventional  CGS  systems  (where  volts  are 
measured  in  the  practical  or  MKS  system).  There  also 
exist  "large"  ions  with  mobilities  between  3 and  80  x 10'4 
cm2/V-s.  These  large  ions  are  presumably  charged  versions 
of  the  so-called  "Aitken  condensation  nuclei.  " The  large 
ions  always  result  from  charge  transfer  collisions  from 
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small  ions.  There  is  also  a class  of  Intermediate  ions  that 
have  been  measured  by  some  observers  under  some  con- 
ditions. However,  it  is  the  Aitken  condensation  nuclei, 
which  exist  in  large  numbers  over  continental  areas  and 
which  have  meteorological  significance,  that  produce  the 
major  complications  in  the  discussion  of  ionic  concentration. 

If  the  effects  of  the  divergence  of  E are  ignored,  then 
the  time  variation  of  the  concentration  of  small  and  large 
ions  can  be  written  as  follows: 


Let:  nj 

"2 

N1 

N2 

N0 

q 

a 

’>12 

>121 

*>10 

”20 

y 


number  concentration  of  positive  small  ions, 
number  concentration  of  negative  small  ions, 
number  concentration  of  positive  large  ions, 
number  concentration  of  negative  large  ions, 

number  concentration  of  Aitken  uncharged 
nuclei, 

ionization  rate  per  unit  volume  for  small  ions, 
recombination  coefficient  for  small  ions, 
recombination  coefficient  for  nj  and  N2, 
recombination  coefficient  for  n2  and  Nj, 
recombination  coefficient  for  nj  and  NQ, 
recombination  coefficient  for  and  N^,  and 
recombination  coefficient  for  Nj  and  N2. 


Then  the  rate  equation  for  positive  small  ions  is 


dnj 

dr  = q-°tnln2'ni2nl  N2  " ”l0  nlN0 


and  a similar  equation  exists  for  n2,  the  negative  small  ion. 
For  the  positive  large  ion: 


dN. 

"dT  = ”10  n10  N0  ‘ ”21  "2  N1  ' yNl  N2  • 
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with  a similar  equation  for  N2.  The  difficulties  of  solving 
this  set  of  equations  are  discussed  in  Ref.  2 and  3.  As  is 
the  isual  case  in  plasma  physics,  the  uncertainties  in  the 
recombination  coefficients  are  so  targe  that  any  results 
are  of  questionable  value.  If,  for  an  equilibrium  solution, 
the  time  derivatives  are  set  equal  to  zero,  and  if  yNj  N2 
and  anj  n2  are  negligible  and  Nj  = N2  = N,  then  it  can  be 
shown  that 


In  other  words,  the  ratio  of  the  concentration  of  Aitken  neu- 
tral ions  to  charged  large  ions  and  the  ratio  of  positive 
small  ions  to  negative  small  ions  depend  upon  the  recom- 
bination coefficients  between  the  small  ions  and  the  charged 
or  uncharged  Aitken  nuclei. 

In  general,  measurements  show  nj/n2  > 1 (1. 1 to  1.  5 
are  typical  values  over  land}  and  nj  varies  from  200  to  600 
cm~3,  the  lower  value  associated  with  high  values  oi  Nq  and 
high  pollution.  Concentrations  of  large  ions  vary  immensely, 
from  200  cm"®  over  the  ocean  to  80  000  cm"®  in  industrial 
areas.  From  Eq.  (9),  using  accepted  rates  of  recombination 
coefficients,  the  value  of  Nq/N  ranges  from  2 to  3.  Under 
any  conditions,  the  atmospheric  conduction  current  is  car- 
ried by  the  small  ions;  the  effect  of  the  Aitken  condensation 
nuclei  is  that  they  affect  the  concentration  of  small  ions. 
These  Aitken  nuclei  are  very  large,  with  radii  of  the  order 
of  2 x 10~6  cm.  Whatever  their  nature,  it  is  well  established 
that  they  are  most  plentiful  in  industrial  areas.  The  life- 
time of  the  small  ion,  as  deduced  from  Eq.  (7),  is  from  20 
seconds  in  industrial  areas  to  300  seconds  over  the  ocean. 
This  large  difference  in  the  lifetime  of  the  small  ion  is  not 
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reflected  in  the  number  concentration  or  conductivity  be- 
cause the  generation  rate  of  small  ions  is  much  greater 
over  land. 

The  above  model  for  the  atmospheric  electric  field 
presents  the  picture  of  a spherical  capacitor,  with  the  outer 
shell  (the  electrosphere)  at  about  +300  000  volts  and  a cur- 
rent flowing  to  the  earth  of  the  order  of  1 x 10"12  A/m2. 
Since  the  typical  value  of  the  potential  gradient  near  the 
earth  is  120  V/m,  from 


a - - t F 
o 


the  charge  density  on  the  surface  of  the  earth  is 
a = -1.C6  x lo'9  C/m2  . 

For  a current  to  the  earth  of  1 x lo*12  A/m2,  the  time  to 
neutralize  the  surface  charge  of  the  earth  is 


1.06  x 10 


= 1060  seconds  or  18  minutes. 


Any  method  of  calculating  the  time  constant  of  atmospheric 
electricity  will  give  a similar  value,  i.  e.,  the  electro- 
sphere-earth condenser  system  will  discharge  in  one-half 
hour  or  less.  If  the  potential  of  the  electrosphere  relative 
to  the  earth  is  to  be  maintained,  some  mechanism  must  be 
found  to  maintain  it.  The  theory  generally  accepted  is  that 
the  atmospheric  potential  is  caused  by  storm  or  thunderstorm 
activity.  Although  the  charging-discharging  balance  has 
never  been  accurately  demonstrated,  it  seems  clear  that 
the  required  flow  of  positive  charge  from  the  earth  occurs 
during  periods  of  electrical  storm  activity.  This  flow  could 
be  due  to  lightning  flashes  or  point  discharges  from  the 
earth  into  the  intense  field-reversed  atmospheric  regions 
that  accompany  severe  electrical  disturbances. 


TMC  JOHN*  HOHM  UMtVMtTT 

APPLIED  PHYSICS  LABORATORY 

•um  , r . m HMnjM 


In  summary,  the  first  order  attempt  at  a model  for 
the  atmospheric  electric  field  posits  a spherical  capacitor 
in  which  the  potential  difference  between  the  electrosphere 
and  the  earth  is  maintained  by  thunderstorm  pumping  and  is 
run  down  by  fair  weather  leakage  currents  elsewhere  on 
the  surface  of  the  earth.  The  serious  deficiencies  in  this 
model  are  well  known  (see  Ref.  4),  but  the  model  remains 
very  useful  for  the  purposes  of  this  investigation. 

The  primary  objection  to  possible  applications  of 
the  atmospheiic  electric  field  is  the  existence  of  fluctua- 
tions of  considerable  magnitude.  Not  only  are  field  rever- 
sals observed  during  times  of  thunderstorm  activity  - and 
sometimes  during  periods  of  precipitation  — but  wide  vari- 
ations in  the  fair  weather  field  are  also  observed.  These 
variations  occur  on  several  time  scales.  There  are  yearly 
and  daily  variations  of  a nonrandom  nature  and  random 
variations  with  time  scales  as  short  as  seconds.  These 
time  variations  show  different  characteristics  over  different 
parts  of  the  earth's  surface.  The  easiest  atmospheric  pa- 
rameter to  measure  is  the  atmospheric  potential  or  electric 
field,  and  most  effort  has  been  expended  on  this  parameter. 

Equation  (3)  can  be  written  as 

F » »L  ' 

where  Rl  is  the  atmospheric  resistivity  at  the  point  of 
measurement.  Since  I = VfR,  where  V is  the  potential  of 
the  electrosphere  and  R is  the  columnar  resistance  of  a 
1 m2  column  of  air  between  the  electrosphere  and  the  earth, 
ien 
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Differentiating  Eq.  (12)  with  respect  to  time  and  dividing 
the  result  by  the  same  equation,  the  result  is 

1 fE.  - i.dV  J_  1_  dR  .... 

F dt  ' V dt  R,  dt  " R dt  ' ' 

Worldwide,  the  yearly  fluctuation  in  the  potential 
gradient  is  a single  cycle,  showing  a maximum  in  January 
and  a minimum  in  July.  Daily,  there  is  a single  cycle 
variation  with  a maximum  at  19  hours  and  a minimum  at  4 
hours,  Greenwich  Mean  Time.  Since  the  correlation  of 
world  thunderstorm  activity  and  potential  gradient  is  very 
good  for  the  daily  potential  gradient  variation,  it  is  reason- 
able to  assign  the  yearly  and  daily  fluctuations  of  the  poten- 
tial gradient  to  the  first  term  in  Eq.  (13),  i.  e. , a variation 
in  the  potential  of  the  electrosphere.  However,  these  sim- 
ple, single-cycle  fluctuations,  which  are  believed  to  be 
characteristics  of  the  average  worldwide  field,  are  only 
observed  over  the  oceans,  the  polar  regions,  and  in  isolated 
continental  regions.  In  the  industrialized  regions  of  the 
world,  where  most  fie'd  measurements  have  been  made 
(e.  g. , England  and  Germany),  the  picture  is  much  more 
complicated.  In  these  locations,  a two-cycle  variation  per 
day  is  often  observed,  which  correlates  with  local  time. 
Thus  atmospheric  pollution,  the  second  term  of  Eq.  (13), 
is  involved.  Term  three  of  Eq.  (13)  is  in-olved  when  R is 
increased  by  a cloud  layer  of  large  vertical  and  horizontal 
extension  that  reduces  the  current  density,  and  consequently 
the  potential  gradient,  if  Rl  remains  a constant. 

In  general,  the  annual  variation  of  V is  about  *15% 
and  daily  variation  about  ±20%.  These  fluctuations,  while 
by  no  means  negligible,  are  small  compared  with  the  fluc- 
tuations caused  by  weather  or  pollution  over  industrialized 
ccntinental  land  areas.  A discussion  of  these  fluctuations, 
which  in  general  are  aperiodic,  is  given  in  great  detail  in 
Chapter  V of  Ref.  3.  Since  the  time  constant  of  the  atmo- 
spheric field  is  about  30  minutes,  the  existence  of  fluctua- 
tions with  periods  of  a few  minutes  or  less  must  be  due  to 
changes  in  the  atmospheric  conditions  near  the  point  of 
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or,  in  terms  of  the  potential  <p, 

V 2<p  + ^ (in  X)  = 0 . 

0 Z 0 z 


If  the  conductivity  is  assumed  to  be  a simple  exponential 
function  of  the  altitude  z. 


X = Xo^z,  then^(tnX)  = 0 . 


and  Eq.  (18)  becomes 

V2V>+  0 = 0 . (19) 

o z 

Two-dimensional  computer  solutions  of  Eq.  (19) 
have  been  made  by  Hoppel  (Eef.  5)  for  a triangular  moun- 
tain and  a cliff.  The  computer  calculations  for  several 
orographic  protrusions  have  been  performed  by  R.  Weiss 
of  API . These  calculations  include  two-dimensional  cal- 
culations for  walls  and  a three-dimensional  calculation  for 
an  upright  pole.  Kwang  Yu  of  APL  has  achieved  an  ana- 
lytic solution  of  Eq.  (18)  for  a class  of  mountain  contours 
(which  happen  to  coincide  very  closely  to  the  contours  of 
typical  mountains  in  the  eastern  U.  S. ).  See  Appendices 
A and  B of  this  report.  Both  Ref.  5 and  Appendix  B con- 
sider effects  caused  by  abrupt  changes  in  the  value  of  the 
atmospheric  conductivity.  All  of  the  above  calculations 
assume  that  the  surface  of  the  earth  and  the  protrusion 
are  perfect  conductors. 

The  equipotential  lines  about  a wall  15  meters  high 
and  2 meters  thick  (computed  by  R.  Weiss)  are  shown  on 
Fig.  1.  The  ground-level  field  value  at  infinity  is  100  V/m 
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and  the  computation  is  two-dimensional,  i.  e. , the  length 
of  the  wall  is  infinite.  While  the  general  shape  of  the  equi- 
potential  surfaces  is  what  one  would  expect,  several  nu- 
merical results  of  the  plots  are  worth  noting.  If  one  comes 
toward  the  wall  along  a plane  5 meters  above  the  ground, 
the  vertical  potential  gradient  decreases  from  80  V/m  30 
meters  from  the  center  of  the  wall  to  19  V/m  4 meters 
from  the  wall.  Along  the  same  plane,  the  horizontal  po- 
tential gradient  is  4 V/m  at  a distance  of  30  meters  and 
30  V/m  at  4 meters.  Moving  in  toward  the  wall  on  a plane 
10  meters  above  the  ground,  the  vertical  potential  gradient 
decreases  from  79  V/m  30  meters  from  the  wall  to  60  V/m 
6 meters  from  the  center  of  the  wall.  Along  the  10-meter 
plane,  the  horizontal  potential  gradient  changes  from  7.  2 
V/m  at  30  meters  to  46  V/m  6 meters  from  the  wall. 
Directly  above  the  wall,  the  vertical  potential  gradient  is 
very  large:  at  a height  of  17  meters  from  the  ground,  the 
vertical  potential  gradient  is  in  excess  of  250  V/m.  The 
calculated  data  show  that  at  a distance  of  two  scale  heights 
from  the  wall,  the  horizontal  gradient  is  rather  small  (less 
than  10  V/m)  and  its  value  is  altitude  dependent.  Obviously, 
close  to  the  ground,  the  horizontal  gradient  approaches 
zero.  However,  two  useful  pieces  of  information  are  avail- 
able as  one  approaches  the  obstacle:  the  horizontal  gradi- 
ent increases  and  the  vertical  gradient  decreases  up  to 
a point,  depending  upon  altitude,  very  close  to  the  wall. 

Weiss'  three-dimensional  calculation  for  a con- 
ducting pole  5 meters  high  and  0. 025  meter  in  diameter  is 
shown  in  Fig.  2.  Obviously,  such  an  orographic  protrusion 
is  very  difficult  to  detect  by  measurements  of  the  perturbed 
atmospheric  field.  On  the  other  hand,  from  the  point  of 
view  of  performing  experiments,  such  a pole  can  be  used 
to  support  instrumentation  without  greatly  perturbing  the 
local  atmospheric  field. 

K.  Yu's  analytic  solution  for  the  equipotential  sur- 
faces in  the  neighborhood  of  South  Mountain,  Maryland  is 
shown  in  Fig.  3 and  discussed  in  Appendix  A.  Because  the 
slope  off  this  mountain  is  typical  of  the  mountains  in  the 
eastern  U.S. , it  is  convenient  for  the  experimental  mea- 
surements discussed  later  in  this  report.  This  calculation 
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was  also  done  for  a nominal  electric  field  value  of  100  V/m. 
At  a distance  of  200  meters  from  the  ridge  at  an  altitude  of 
350  meters,  the  horizontal  field  is  33  V/m.  One  thousand 
meters  from  the  ridge  at  an  altitude  of  275  meters,  the 
horizontal  field  is  5.  5 V/m.  Figure  A4  of  Appendix  A 
shows  the  results  of  a constant  contour  plot  for  the  ampli- 
tude of  the  electric  field.  If  the  mountain  is  approached  on 
a horizontal  plane  lying  below  the  mountain  peak,  the  elec- 
tric field  amplitude  |£  | decreases  until  at  some  point  (de- 
pending upon  altitude)  the  field  begins  to  increase.  At  an 
altitude  of  300  meters  the  field  minimum  occurs  about  450 
meters  from  the  ridge. 

The  calculations  for  South  Mountain  show  that  at  a 
horizontal  range  of  1000  meters  the  horizontal  gradient  mea- 
sured at  275  meters  above  sea  level  (110  meters  above  the 
ground)  is  5.  5 V/m.  This  mountain  is  375  meters  above 
sea  level.  The  height  of  the  mountain,  measured  relative 
to  ground  level  at  the  range  of  1000  meters  from  the  peak, 
is  210  meters.  Thus  the  horizontal  gradient  at  the  range 
of  five  times  the  height  from  the  mountain  is  computed  to 
be  5. 5 V/m  for  an  altitude  of  about  one-half  the  mountain 
height.  This  result  is  sensitive  to  the  contour  of  the  pro- 
trusion; for  the  wall  of  Fig.  2,  the  horizontal  gradient  is  at 
most  1. 5 V/m  at  a range  of  five  times  its  height  from  the 
wall.  Further,  because  of  the  second  term  of  Eq.  (19),  the 
solution  of  Eq.  (19)  for  any  protrusion  must  be  scaled  with 
caution.  Since  0 “■  2.  5 * 10'*  m'1,  the  change  in  conduc- 
tivity with  altitude  is  of  little  significance  for  low  protru- 
sions. Therefore,  if  for  the  same  field  strength,  the  15- 
meter  wall  of  Fig.  2 is  replaced  by  a wall  5 meters  high, 
the  vertical  and  horizontal  gradients  will  be  essentially 
identical  if  measured  at  ranges  and  heights  divided  by  a 
factor  of  3.  For  the  mountain  in  Fig.  3,  if  the  height  is 
increased,  say  by  a factor  of  2,  and  the  geometric  similar- 
ity of  the  slope  is  preserved,  then  under  the  assumptions 
made  in  Appendix  A,  the  shape  of  the  equipotential  surfaces 
is  preserved.  However,  value  in  volts  of  these  equipoten- 
tial surfaces  will  not  precisely  double  for  a point  of  geo- 
metric similarity  if  the  unperturbed  field  strength  remains 
the  same.  While  the  effects  of  scale  changes  on  the  values 
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of  the  vertical  and  horizontal  gradients  have  not  been  nu- 
merically computed,  they  will  certainly  be  small  over  the 
altitude  range  for  which  Eq.  (5)  is  valid.  This  range  in- 
cludes all  typical  mountains  in  the  eastern  U.S. 

In  the  above  calculations  an  exponential  conductivity 
of  the  form  X = X0e^z  is  used.  A question  arises  as  to  the 
validity  of  this  assumption  for  orographic  protrusions  such 
as  mountains.  Hoppel  (Ref.  5)  has  considered  the  effect  of 
the  change  in  conductivity  that  occurs  at  an  inversion  layer 
at  1400  or  3500  meters  on  his  triangular  mountain  calcula- 
tion for  a mountain  3400  meters  high  by  using  different 
values  of  XQ  and  p above  and  below  the  inversion  layer. 

On  the  realistic  assumption  that  \0  and  lip  increase  above 
the  inversion  layer,  the  effect  is  to  increase  the  curvature 
of  the  equipotential  lines  above  the  protrusion,  which  would 
tend  to  aid  the  detection  of  the  presence  of  the  mountain. 


In  Appendix  B,  K.  Yu  considers  the  problem  of  a 
local  semicircular  region  of  higher  conductivity  located  on 
the  top  of  a mountain  ridge.  This  would  correspond  to  a 
situation  when  ground  radioactivity  increased  the  conduc- 
tivity at  the  top  of  the  ridge  or  where  pollution  or  humidity 
levels  decreased  the  conductivity  in  the  valley.  The  result 
is  a sharp  jump  upward  in  the  equipotential  lines  at  the 
point  of  discontinuity  in  the  conductivity.  While  no  such 
sharp  changes  in  conductivity  could  be  expected  in  nature, 
the  conclusions  are  qualitatively  correct.  One  would  expect 
that  the  conductivity  in  the  valleys  (if  it  departs  from  the 
idealized  exponential  behavior)  will  be  less  than  expected 
because  of  pollutants  and  thus  compress  the  equipotential 
lines.  In  general,  the  slope  of  the  equipotential  lines  lead- 
ing into  the  mountain  will  be  increased  and  the  horizontal 
gradient  will  be  increased. 

In  summary  (from  the  point  of  terrain  avoidance), 
it  is  reasonable  to  expect  that  the  departure  of  the  atmo- 
spheric conductivity  from  the  simple  exponential  law  can 
make  the  detection  of  the  horizontal  gradient  easier. 
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3.  METHODS  OF  MEASURING  THE 
ATMOSPHERIC  ELECTRIC  FIELD 


The  moe  t convenient  atmospheric  electric  field 
parameter  to  measure  is  either  the  atmospheric  potential 
or  the  electric  Held.  While  the  atmospheric  electric  field 
is  of  the  order  of  100  V/m  near  the  ground,  the  problem  of 
measuring  it  is  made  very  difficult  by  the  fact  that  the  at- 
mospheric resistivity  is  typically  5 x lo13  ohms/m.  Thus 
the  atmospheric  potential  (relative  to  some  known  reference 
which  may  be  earth)  requires  a very  high  impedance  volt- 
meter, or  a true  electrostatic  voltmeter.  The  electric 
field  may  be  measured  directly  using  the  relation 

E - O l(  , (20) 

o 

where  o is  the  bound  charge  density  on  a conductor  in  the 
electric  field.  In  practice,  it  is  necessary  to  mechanically 
modulate  the  field  on  the  conductor  and  measure  the  transfer 
of  charge  (i.  e. , current)  during  the  modulation  process. 

Such  a device  is  called  a field  mill. 

The  history  and  methods  of  measuring  the  atmo- 
spheric electric  potential  or  field  are  reviewed  in  Refs.  2 
and  3.  The  development  of  modern  solid  state  electronics 
has  ameliorated,  if  not  eliminated,  the  difficulties  of  mea- 
surement. For  these  experiments,  two  devices  were  used, 
a high  impedance  voltmeter  coupled  to  the  atmosphere  with 
a radioactive  equalizer  and  a cylindrical  field  mill. 

The  high  impedance  voltmeter  used  in  this  exper- 
iment was  designed  at  APL  by  R.  Konigsberg  (Ref.  6). 


Ref.  6.  R.  L.  Konigsberg,  "Design  of  Portable, 
Solid  State,  High  Input  Impedance  High  Voltage  Differential 
Electrostatic  Voltmeter,  " submitted  to  Electronic  Design, 
13  January  1975. 
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This  instrument  uses  two  radioactive  atmospheric  couplers 
or  sensors,  usually  placed  a meter  apart.  If  Vj  and  V2 
are  the  voltages  measured  at  the  sensor  positions,  the 
signals  are  combined  to  yield: 

Vj  - Vg  = differential  voltage,  0 to  500  volts 

Vl+ V2 

= common  mode  voltage,  0 to  250  volts. 

Figure  4 is  a block  diagram  of  this  instrument  All  output 
voltages  are  referenced  to  instrument-case  ground,  and 
the  dynamic  range  of  the  instrument  is  limited  by  the  high 
voltage  supply  and  die  high  voltage  emitter  follower.  The 
input  impedance  of  the  voltmeter  is  10*4  ohms  or  greater. 
The  input  capacity  is  reduced  by  feeding  the  output  of  the 
emitter  follower  (L  e. , the  sensor  potential)  back  to  the 
cable  shielding  the  input  and  the  guard  plate  below  the  radio- 
active equalizer.  For  the  two  instruments  that  have  been 
constructed,  the  input  capacity  measures  about  1 pF.  The 
details  of  the  radioactive  sensor,  guard  circuit,  and  mount- 
ing arm  are  shown  on  Fig.  5. 

The  radioactive  sensor  used  on  this  instrument  has 
usually  been  a 500  pCi  polonium  210  Ct  emitter.  Since  the 
range  of  these  a particles  in  air  is  about  4 cm,  a relatively 
intense  region  of  ionization  surrounds  the  potential  equal  - 
izer.  The  effective  couplirg  impedance  of  such  equalizers 
has  never  been  precisely  calculated,  although  values  of 
2 x 1010  ohms  have  been  quoted  (Ref.  2).  For  the  present 
instrument,  if  the  potential  at  the  point  read  by  the  radio- 
active instrument  is  to  be  accurate,  the  coupling  impedance 
to  the  atmosphere  obviously  must  be  much  less  than  10l4 
ohms.  Since  the  time  response  of  this  instrument  in  the 
atmospheric  field  or  a laboratory  Reid  is  measured  to  be 
less  than  1/2  second  (and  is  probably  faster  than  the  re- 
corder used  in  the  measurement)  from  the  1 pF  input  ca- 
pacity, the  coupling  impedance  is  certainly  less  than  5 x 101 1 
ohms.  Therefore,  the  inherent  error  of  this  instrument  in 
measuring  the  atmospheric  potential  is  negligible. 
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The  cylindrical  field  mill  used  in  this  experiment 
was  one  constructed  by  H.  W.  Kasemir  (Ref.  7)  and  loaned 
to  APL  by  the  Army  Electronics  Command  (ECOM),  Ft. 
Monmouth.  While  simple  in  principle,  the  field  mill  is 
complicated  mechanically  and  expensive  to  construct.  Fig- 
ure 6 is  a greatly  simplified  schematic  of  the  cylindrical 
field  mill.  The  mill  consists  of  an  insulated  split  cylinder. 
This  cylinder  is  supported  by  a base  and  rotated  by  a small 
motor  within  the  base.  Using  an  infinite  cylinder  approxi- 
mation, the  charge  on  the  top  segment  of  the  cylinder  is 
given  by  (Ref.  6) 

Q = -4«  EaL  sin  V.  , 
o l 

where  E is  the  electric  Reid,  a and  L the  radius  and  length 
of  the  cylinder,  and  <Pi  the  angle  between  the  vector  fi  and 
the  insulated  strip  between  the  segments  of  the  cylinder. 

If  the  cylinder  rotates  with  an  angular  velocity  a:,  then  the 
current  flowing  between  the  segments  is 

I = -4<EaLut  cos  (n-t  + (Pj)  . 

The  voltage  output  of  an  operational  amplifier  with  a trans- 
fer impedance  of  R,  mounted  between  the  cylindrical  seg- 
ments is 

U = -IR  = 4fEaLo;R  cos  tet  + (p.)  . (21) 

o 1 

For  this  field  mill,  a = 5 cm,  L = 10  cm,  R = 22xlo6 
ohms,  and  the  frequency  of  revolution  is  1800  rpm.  Thus 
for  an  electric  field  of  100  V/m  the  peak  value  of  UQ  is  77 
mV.  The  output  of  the  operational  amplifier  (mounted  in- 
side the  revolving  cylinder)  is  taken  from  a slip  ring  and 
rectified  by  a phase  sensitive  rectifier  locked  to  the  angle 
4>1.  The  rectified  signal  is  amplified  and  recorded.  A 
complicated  set  of  slip  rings  (not  shown  on  Fig.  S)  is 


Ref.  7.  H.  W.  Kasemir,  The  Cylindrical  Field 
Mill,  AD-610485,  ECOM,  Ft.  Monmouth,  N.  J.,  1964. 
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required  to  supply  power  to  the  current  amplifier  and  pro- 
vides the  phase  information  for  the  rectifier.  By  using  two 
phase-locked  amplifiers,  phased  90“  apart,  both  the  x and 
z components  of  the  electric  field  can  be  recorded  simul- 
taneously. The  calibration  of  the  field  mill  and  the  orien- 
tation of  its  base  to  give  orothogonal  vertical  and  horizon- 
tal field  measurements  were  performed  in  the  laboratory 
by  placing  the  mill  between  1 m2  capacitor  plates  1.64 
meters  apart.  The  dynamic  range  of  the  field  mill  is  very 
large:  for  a single  value  of  R * 22  megohms,  the  field  mill 
will  record  field  values  of  from  25  to  1250  V/m  full  scale, 
and  the  value  of  R can  be  remotely  switched  to  record  fields 
of  the  order  of  100  000  V/m.  Thus,  with  a good  recorder, 
field  values  as  low  as  0.  25  V/m  can  be  recorded.  Unfor- 
tunately, the  z and  x gain  channels  of  this  field  mill  cannot 
be  set  at  separate  values  with  the  present  electronics.  How- 
ever, overdriving  one  channel  does  not  affect  the  other. 


Whatever  method  is  used  to  measure  the  atmo- 
spheric field,  the  instrument  must  be  supported  in  the  at- 
mosphere and  the  data  recorded.  This  support  system 
necessarily  produces  some  distortion  in  the  atmospheric 
field.  Further,  the  observer  must  remain  some  distance 
from  the  measuring  apparatus  since  he  is  an  orographic 
protrusion  and  his  clothing  is  usually  electrostatically 
charged. 


For  measurements  with  the  high  impedance  volt- 
meter, the  support  stand  shown  in  Fig.  7 was  used.  The 
vertical  support  rods  were  always  grounded;  the  horizontal 
bar  supporting  the  voltmeter  was  either  grounded  or  op- 
erated at  a positive  bias  with  respect  to  ground.  If  the 
positive  bias  is  set  to  a value  so  that  the  common  mode 
voltage  of  the  meter  is  zero,  then  the  support  bar  and  meter 
case  are  at  local  atmospheric  potential  and  very  little  per- 
turbation of  the  local  field  occurs.  Meter  readings  were 
either  taken  from  a distance  with  a telescope  or  recorded 
on  an  accurate  potentiometer-type  recorder  mounted  on 
the  ground  15  or  20  feet  away  from  the  instrument.  When 
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Fig.  8 Position  of  Atmospheric  Potential  Field  Sensors  on  Cessna  337  Airplane 
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the  case  of  the  instrument  is  biased  to  atmospheric  poten- 
tial, the  output  signal  of  the  meter  must  be  isolated 
from  ground.  The  recorder  used  for  these  measurements 
has  a double-ended  input  that  provides  about  600  volts  iso- 
lation from  ground.  The  recorder  accurately  indicates  the 
difference  voltage  applied  across  its  terminal.  A signal 
isolator  was  constructed  and  used  occasionally  for  con- 
verting the  DC  meter  output  into  a pulse-width  modulated 
signal.  This  signal  was  used  to  drive  a light-emitting  diode 
diode.  The  light  signal  was  detected  and  integrated,  yield- 
ing a signal  referenced  to  ground  that  was  identical  to  the 
input  signal.  This  device  provided  isolation  for  at  lea&i 
1500  volts.  In  general,  with  a ground-based  recorder,  one 
is  limited  to  heights  of  less  than  10  meters. 

Aircraft  measurements  of  the  atmospheric  field 
were  made  with  an  instrumented  Cessna  337  aircraft.  Ra- 
dioactive potential  equilizers  were  mounted  on  the  wing  tips 
and  tail  rudders  (Fig.  8).  Each  sensor  was  mounted  about 
0. 5 meter  from  the  metal  skin  of  the  aircraft.  Sensors  1 
through  6 were  connected  to  the  high  impedance  voltmeter 
circuit  (Fig.  4),  and  the  signals  could  be  combined  in  the 
following  modes:  (Vl  + V2>/2,  Vi  - V2,  (V3  J-  V4)/2, 

V3  - V4,  (VX  + V3)/2,  Vi  - V3,  (V2  + V4)/2,  V2  - v4, 

(V5  + Vg)/2,  and  Vg  - Vg.  Sensors  7 and  8 were  connected 
to  a low  impedance  autopilot  circuit  (described  in  Ref.  1) 
and  could  be  combined  to  yield  (V7  + Vg)/2  and  V7  - Vg. 
During  flight  the  processed  sensor  signals  were  recorded 
on  magnetic  tape  (20  channels  available),  and  two  channels 
were  available  on  a strip  recorder.  For  the  flights  cov- 
ered in  this  report,  not  all  of  the  above  signal  combinations 
were  available  at  one  time  because  of  a lack  of  operational 
amplifier  summing  and  differencing  channels.  For  example, 
for  the  wing  sensors  (Vi  + V2>/2,  Vi  - V2,  (V3  + V4)/2, 
and  V3  - V4  could  be  recorded  simultaneously.  By  throwing 
a switch,  the  signal  processing  could  be  changed  to  record 
(Vi  + V3)/2,  Vi  - V3,  <V2  + V4)/2,  and  V2  + V4.  The  air- 
craft was  also  equipped  with  a vertical  gyro  so  that  bank 
and  pitch  angles  could  be  recorded.  The  pods,  which  con- 
tained instrumentation  to  measure  conductivity,  were 
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mounted  under  each  wing.  Various  flight  parameters  such 
as  air  speed,  altitude,  and  temperature  could  also  be  re- 
corded, However,  the  aircraft  instrumentation  was  pri- 
marily dedicated  to  the  ECOM -sponsored  RPV  program; 
for  the  purposes  of  terrain  avoidance,  the  signals  of  pri- 
mary importance  were  from  the  wing  and  tail  atmospheric 
potential  sensors. 

Measurement  of  the  atmospheric  potential  from  an 
aircraft  faces  two  difficulties:  (a)  the  potential  at  the  sen- 
sor position  must  be  referenced  to  something,  usually  the 
aircraft  skin,  and  (b)  there  is  a distortion  in  the  potential 
field  caused  by  the  presence  of  the  conducting  skin  of  the 
aircraft. 

It  is  well  known  that  the  aircraft  skin  can,  under 
certain  conditions,  become  charged  to  very  high  potentials 
(hundreds  of  thousands  of  volts)  relative  to  the  ambient 
atmospheric  potential.  The  perturbation  of  the  potential 
field  caused  by  a conductor  can  be  calculated  analytically, 
but  these  calculations  are  limited  to  those  few  simple 
shapes  for  which  Laplace's  equation  can  be  solved.  As  a 
consequence,  the  "form  factors"  for  the  potential  field  or 
the  "augmentation  factors"  for  the  electric  field  must  be 
determined  experimentally. 


The  difficulties  of  measuring  the  atmospheric  po- 
tential from  an  aircraft  are  closely  related  to  the  problem 
of  the  electrostatic  autopilot.  This  situation  is  covered  in 
Appendix  C,  where  the  potentials  and  fields  in  the  neigh- 
borhood of  conducting  spheres  and  prolate  spheroids  are 
discussed.  For  the  case  of  the  conducting  sphere,  using 
the  coordinate  system  of  Appendix  C,  the  potential  outside 
the  sphere  is  given  bi  * 
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I 
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where  <p0  is  the  atmospheric  potential  and  <PS  is  the  poten- 
tial on  the  surface  of  the  sphere  caused  by  self-charge.  If 
a voltage  sensor  is  placed  outside  the  sphere,  reading  the 
atmospheric  potential  relative  to  the  skin  of  the  sphere, 
then  that  sensor  will  read  a voltage  given  by 


P ft  - \ ) + P <%  - 1)  • 

o J s r 
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The  case  of  the  prolate  spheroid  leads  to  results  of  a simi- 
lar form  and,  in  general,  a sensor  mounted  off  the  skin  of 
the  aircraft  will  measure  a voltage  given  by 


*01  <ri>  - D 


(22) 


where  (?},  which  is  a function  of  the  aircraft  surface 
shape  and  the  location  of  the  sensor,  is  the  perturbation 
caused  by  the  charge  induced  on  the  aircraft  skin  by  the 
atmospheric  field;  Pi  (?i)  is  the  similar  factor  for  the  po- 
tential field  caus  <d  by  any  self-charge  on  the  aircraft. 

Pi  (ri)  is  a fixed  function  for  any  coordinate  system  fixed 
in  the  aircraft,  but  Ot i (?)  is  a function  of  aircraft  position. 
Calculations  for  the  sphere  show  (if  the  analytic  expression 
for  the  potential  is  derived  for  <PQ  - 0 at  the  center  of  the 
sphere)  that  convention  for  the  assigning  of  values  for  the 
atmospheric  potential  must  be  preserved  for  all  calculations. 

It  is  clear  from  Eq.  (22)  that  for  a high  impedance 
voltmeter  coupled  to  the  atmosphere  with  a radioactive 
potential  equalizer  (that  reads  the  atmospheric  potential 
at  its  location  relative  to  <Pg),  if  Ps  is  large  enough  the  sig- 
nal Vj  will  saturate  the  amglifier.  (For  sensors  located 
off  the  skin,  oij  (r)  and  0 j (r)  are  always  less  ihan  one.) 
Furthermore,  because  the  dynamic  range  of  this  high  im- 
pedance instrumentation  is  V « 250  volts  for  any  sensor, 
the  atmospheric  field  is  often  greater  than  100  V/m  and 
the  dimensions  of  the  aircraft  are  several  meters,  the 
point  in  the  atmosphere  at  which  Ps  coincides  with  the  at- 
mospheric potential  is  very  important.  For  example,  if 
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the  atmospheric  potential  coincides  with  the  skin  potential 
for  a point  low  on  the  aircraft  skin,  the  sensors  mounted 
above  the  top  of  the  wings  can  easily  saturate  their  ampli- 
fiers. If  Vj  saturates,  no  signal  processing  such  as 
Vi  - V2  can  yield  a valid  result,  even  though  V]  - V2  is 
well  within  the  dynamic  range  of  the  instrument. 


Early  in  the  history  of  the  flights  with  the  Cessna 
337,  it  was  realized  that  such  saturation  effects  did  indeed 
occur.  The  analysis  of  experimental  data  was  consistent 
with  the  fact  that  the  aircraft  skin  potential  assumed  an 
atmospheric  potential  for  a point  low  on  the  aircraft  frame. 
Presumably  this  point  coincides  closely  with  the  engine 
exhausts.  In  general,  for  fair  weather  flights,  no  high 
values  of  ipB  have  been  observed.  It  has  been  necessary, 
however,  to  provide  a method  of  biasing  the  aircraft  frame 
so  that  signals  such  as  Vj  - V3  could  be  measured  over  a 
range  of  bank  angles.  This  was  accomplished  by  mounting 
a conducting  sheet  to  the  aircraft  skin  behind  the  exhaust 
of  the  front  engine.  For  example,  by  applying  up  to  1000 
volts  on  this  conductor  it  is  possible  to  set  the  skin  po- 
tential so  that  Vj  is  nearly  zero,  or  if  the  experiment  re- 
quires it,  making  Vj  positive  and  V2  negative.  The  only 
way  to  assure  that  saturation  does  not  occur  is  to  measure 
the  output  of  each  sensor  individually  before  it  is  pro- 
cessed. For  each  flight  the  skin  bias  must  be  adjusted  to 
a value  to  assure  nonsaturation  for  all  the  sensors  relevant 
to  the  test  plan  for  that  flight. 

For  most  of  this  experimental  work,  both  on  the 
ground  and  in  the  air,  the  high  impedance  voltmeter  coupled 
to  the  atmosphere  with  a radioactive  potential  equalizer 
has  been  used.  Problems  inherent  in  the  use  of  radioactive 
equalizers  have  been  recognized  for  some  time  (Refs.  2,  3, 
and  8).  Unfortunately,  no  general  theory  of  the  radioactive 
potential  equalizer  exists.  Recently,  M.  L.  Hill  and 

Ref.  8 R.  Wagner,  "Zur  Messung  des  luftelekrishen 
Potentialgefalles  mittels  Kollektoren,  " Archiv  fur 
Meteorologic,  Geophvsik  und  Bioklimatologic,  Series  A, 

Vol.  3,  No.  4,  1955,  pp.  427-463. 
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W.  A.  Hoppel  (Ref.  9)  of  APL  analyzed  the  performance  of 
the  radioactive  equalizer  to  a greater  degree  than  had  been 
previously  done.  Their  emphasis  was  placed  primarily  on 
the  autopilot  circuit  for  RPVs.  These  are  relatively  low 
impedance  circuits,  and  since  the  coupling  of  the  sensor  to 
the  atmosphere  is  a higher  impedance  than  that  presented 
by  the  amplifier,  the  autopilot  circuit  and  the  radioactive 
equalizer  are  best  treated  in  a current  source  mode.  The 
difference  signal  between  the  right  and  left  wing  sensors, 
for  example,  is  approximately  given  by 

i,  -L,  = 2GAE  sin  8 , 

L R z 


\ 


where  8 is  the  bank  angle  of  the  aircraft,  Ez  is  the  vertical 
electrical  field,  A is  the  electric  field  augmentation  factor, 
and  G is  the  gain  of  the  electric  field  sensors  and  the  as- 
sociated electronics.  The  conclusions,  which  are  well 
supported  by  experiment,  are  that  the  current  flow  through 
the  sensor  is  determined  by  the  atmospheric  electric  field, 
the  radioactivity  of  the  at  source,  the  dimensions  of  the 
source,  and  the  air  velocity  flowing  over  the  source.  The 
theory  does  not  work  for  zero  air  velocity,  and  for  the 
i ealistic  cases  they  consider,  the  current  tends  to  become 
independent  of  air  velocity  for  air  speeds  above  20  mi/h. 
Thus  they  substantiate  what  numerous  RPV  flights  have 
shown,  that  the  equalizer  in  its  current  operating  mode  is 
a reliable  indicator  for  the  atmospheric  field  as  long  as 
the  radioactivity  is  not  reduced  by  water  or  oil  on  the  sur- 
face of  the  at  emitter. 


For  the  case  of  the  radioactive  equalizer  connected 
to  the  high  impedance  voltmeter.  Hill  and  Hoppel  (Ref.  9) 
show  that,  for  a cylindrically  symmetrical  ionizer,  the 


Ref.  9.  M.  I.,  hill  and  W.  A.  Hoppel,  “Effects  of 
Velocity  and  Other  Physical  Variables  on  the  Currents  and 
Potentials  Generated  by  Radioactive  Collectors  in  Electric 
Field  Measurements,  " presented  at  the  Fifth  International 
Conference  on  Atmospheric  Electricity,  Garmisch- 
Partenkirchen,  West  Germany,  2-7  September  1974. 
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ionizer  assumes  the  corree*  atmospheric  potential  provided 
the  net  current  drawn  from  the  sensor  is  essentially  zero. 
This  result  is  independent  of  the  air  velocity  over  the  sen- 
sor. If  the  sensor  is  not  symmetrical  about  its  axis  per- 
pendicular to  the  atmospheric  field,  a positional  error  is 
produced.  This  error  is  of  the  order  of  the  displacement 
produced  by  the  highly  conducting  sheath  about  the  radio- 
active probe.  In  practice  it  is  relatively  small  since  the 
high  conductivity  region  extends  for  only  a few  centimeters 
They  further  show  that  if  the  equalizer  has  an  initial  charge 
or  if  one  is  imparted  to  it,  the  time  constant  for  the  equal- 
izer to  reach  equilibrium  is 


r-  e0/X 


(23) 


where  X is  the  conductivity  in  the  highly  conducting  region. 
Since  in  this  region  X is  orders  of  magnitude  greater  than 
its  atmospheric  value  of  10-14  mho/m,  Tcan  be  expected 
to  be  less  than  1 second.  From  Ref.  9,  a conservative 
estimate  of  r is  0. 2 second  for  a 500  MCi  a source.  T is  a 
function  of  the  radioactivity  of  the  source:  however,  the 
voltmeter  reading  is  not  as  long  as  the  sensor  coupling 
impedance  to  the  atmosphere  is  much  less  than  the  input 
impedance  of  the  meter. 
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4.  RESULTS  OF  THE  EXPERIMENTAL  PROGRAM 


i 


NOISE  ON  THE  HIGH  IMPEDANCE  VOLTMETER 


When  the  high  impedance  voltmeter  was  mounted  on 
the  support  structure  and  the  common  mode  and  differen- 
tial signals  recorded,  one  of  the  surprising  results  was 
the  high  noise  level  on  the  differential  signal  (see  Fig.  7). 
The  voltmeter  was  usually  mounted  with  the  arms  sup- 
porting the  radioactive  sensors  in  a horizontal  position, 
consequently  the  common  mode  signal  ((Vj  + V2>/2)  was  a 
measure  of  the  vertical  field  and  the  differential  signal  mea- 
sured the  horizontal  gradient.  Figure  9 shows  an  example 
of  the  fluctuations  observed.  These  data  were  recorded  on 
the  flat  roof  of  an  APL  building  in  Howard  County,  Maryland; 
the  departure  of  the  average  differential  readings  from  zero 
is  explained  by  the  fact  that  the  voltmeter  was  mounted 
close  to  one  edge  of  the  flat  roof. 


The  Fig.  9 data  were  recorded  on  consecutive  days 
and,  even  though  the  field  strengths  are  comparable,  the 
fluctuations  for  the  second  day  are  much  greater;  and  ex- 
ceed ±10  volts  about  the  average  for  the  differential  chan- 
nel. The  weather  conditions  were  different  and  may  be 
summarized  as  follows 


Temperature  (°C) 

6 June  1974 

7 June  1974 

26 

15.  5 

Relative  humidity  (%) 

30 

69 

Wind  velocity  (mi/h) 

0-2 

4-5,  gusty 

Barametric  pressure  (mb) 

1010 

1011 

Cloud  cover 

Clear,  sunny  day 

Overcast,  threatening 

It  is  difficult,  however,  to  find  any  reliable  correlation 
between  meteorological  data  and  the  measured  fluctuations. 
Figure  10  shows  data  that  were  recorded  on  a clear  day. 
The  wind  velocity  was  from  10  to  20  mi/h.  These  data 
are  vary  noisy.  The  data  shown  in  Fig.  9 were  recorded 
with  the  sensors  1 meter  above  the  roof  surface;  for  the 
data  shown  in  Fig.  10  the  sensors  were  1.  36  meters  high, 
and  thus  the  average  vertical  electrical  field  is  about  37% 
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higher.  Figure  10  is  recorded  on  a faster  time  scale  to 
show  the  high  speed  (<  1 second)  of  the  fluctuations. 

Battery  voltages  applied  to  the  sensors  through  10^ 
ohm  resistors  proved  that  the  battery-powered  high  im- 
pedance voltmeter  has  negligible  internal  noise:  read- 
ings taken  in  open  fields  showed  that  the  above  fluctuations 
are  not  caused  by  the  local  laboratory  environment.  Read- 
ings taken  with  the  sensors  at  altitudes  of  from  1 to  4.  5 
meters  show  no  significant  difference.  The  fluctuations  are 
definitely  coupled  into  the  voltmeter  through  the  radioactive 
sensors,  but  the  question  remained  as  to  whether  the  fluc- 
tuations were  characteristic  of  the  atmospheric  electric 
field  or  were  the  result  of  some  interaction  between  the 
radioactive  potential  equalizer  and  the  atmosphere. 

The  data  shown  in  Figs.  9 and  10  were  taken  with 
the  high  impedance  voltmeter  case  grounded;  that  is,  the 
voltmeter  was  connected  to  the  recorder  by  shielded  cable 
with  the  cable  shield  connected  to  the  instrument  case  and 
grounded  at  the  recorder.  This  distorts  the  atmospheric 
field  about  the  instrument  and  produces  an  error  in  the  at- 
mospheric potential  reading  at  the  sensor  position.  Iso- 
lating the  recorder  input  from  ground  and  biasing  the 
voltmeter  case  to  near  atmospheric  potential  increases  the 
atmospheric  potential  as  read  at  the  sensor  and,  for  all  the 
cases  in  which  this  has  been  done,  decreases  the  differential 
fluctuations.  An  example  of  these  effects  is  shown  in 
Fig.  11.  The  sensors  were  1.  27  meters  above  the  flat  roof 
top.  For  the  grounded  voltmeter,  the  common  mode  reads 
about  105  volts;  when  biased  to  near  atmospheric  potential 
(i.  e. , common  mode  near  zero),  the  instrument  bias  reads 
about  155  volts.  Thus  the  high  impedance  voltmeter,  when 
operated  with  its  case  grounded,  reads  about  67%  of  the  true 
atmospheric  potential.  While  biasing  the  high  impedance 
voltmeter  to  atmospheric  potential  does  not  eliminate  the 
fluctuations  in  the  horizontal  gradient  (the  differential  sig- 
nal), the  reduction  in  their  amplitude  has  always  been  found 
to  be  significant. 
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In  an  attempt  to  determine  the  origin  of  these  fluc- 
tuations, a series  of  laboratory  measurements  was  made  ; j 

with  the  high  impedance  voltmeter.  In  the  laboratory,  the 
atmospheric  electric  field  does  not  exist.  To  generate  a 
field,  the  meter  sensors  were  placed  between  1 m2  metal  1 

plates,  1.  64  meters  apart, and  the  top  plate  was  connected 
to  a precision  high  voltage  supply.  When  voltage  was  ap- 
plied between  the  plates,  fluctuations  appeared  on  the  volt- 
meter outputs  and  increased  in  amplitude  as  the  electric 
field  was  increased.  The  fluctuations  were  similar  to  those  - , 

observed  in  the  atmosphere,  except  that  their  frequency  J 

w»s  somewhat  lower.  In  the  laboratory  the  potential  across 
the  capacitor  plates  does  not  fluctuate  and  the  electric  field  . 

in  the  absence  of  the  sensor  cannot  fluctuate.  It  was  sus-  ; 

pected  that  the  fluctuations  recorded  above  are  (at  least  in 
part)  spurious  in  the  sense  that  they  do  not  measure  the  at- 
mospheric field  accurately.  J 


This  was  confirmed  by  comparing  the  high  impedance 
voltmeter  with  the  cylindrical  field  mill.  Figure  12  shows 
a comparison  of  field  mill  and  the  high  impedance  voltmeter 
readings  when  they  were  placed  between  capacitor  plates  in 
the  laboratory.  While  the  high  impedance  voltmeter  shows 
its  characteristic  fluctuations,  the  field  mill  is  essentially 
free  of  fluctuations.  The  difference  of  scale  shown  on 
Fig.  12  i3  not  significant;  the  field  mill  and  its  support 
were  calibrated  between  the  capacitor  plates  and  the  mea- 
sured electric  field  is  correct  by  definition.  The  high  im-  ■ I . 

pedance  voltmeter  reads  potential,  not  Held,  and  its  reading  i ' < 

depends  on  how  high  it  is  placed  from  the  bottom  capacitor 
plate.  An  additional  uncertainty  is  produced  because  the 
sensors  extend  into  the  fringing  region  of  the  capacitor 
plates.  Thus  it  would  appear  that  the  fluctuations  are 
caused  by  some  interaction  of  the  radioactive  source  and 
the  air  when  in  an  electric  field.  This  is  proven  in  Fig.  13, 
where  the  vertical  and  horizontal  field  strengths  measured 
by  the  field  mill  are  recorded  with  and  without  a 500  ^Ci  oi 
source  suspended  10  cm  from  the  side  of  the  field  mill. 

When  the  field  mill  was  mounted  outside  in  the  atmosphere 
it  did  not  show  the  fluctuation  characteristic  of  the  high  im- 
pedance voltmeter.  Figure  14  illustrates  the  readings  of 
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the  field  mill  taken  on  7 June  1974  and  is  to  be  compared 
with  the  high  impedance  voltmeter  reading  of  the  same  date  ! 

shown  in  Fig.  9.  The  vertical  field  shows  appreciable  *■ 

changes  over  periods  of  several  minutes,  but  there  is  no  1 

measured  horizontal  field.  As  pointed  out  in  Section  3,  i 

the  time  response  of  the  field  mill  to  changes  in  the  electric  ' ■ 

field  is  faster  than  0.  5 second,  and  therefore  the  smooth- 
ness of  the  above  field  mill  recordings  is  not  the  result  of  1 ; 

instrumental  filtering.  1 , 

From  the  c^'scussion  of  the  atmospheric  field  about  . I :• 

orographic  protrusions  in  Section  2,  it  is  clear  that  the  ‘ ‘ j 

fluctuations  observed  with  the  high  impedance  voltmeter  1 

coupled  to  the  atmosphere  with  radioactive  sensors  pose  a 
serious  problem  since  the  fluctuations  are  of  the  same  or- 
der of  magnitude  as  the  potential  gradients  that  must  be  de- 
tected. In  particular,  ground-based  measurements  of  the 
atmospheric  field  about  orographic  protrusions,  which  are 
essentially  static  measurements,  will  be  subject  to  many 
inaccuracies  and  the  application  of  the  high  impedance 
voltmeter  for  terrain  avoidance  to  slowly  moving  vehicles 
will  be  very  difficult.  However,  one  might  expect  that  if 
the  air  stream  moves  over  the  sensor  with  sufficiently  high 
velocity,  the  fluctuations  will  be  reduced  or  eliminated. 

The  radioactive  sensor  is  surrounded  by  a region  of 
relatively  highly  ionized  and  conducting  air.  The  position 
of  this  highly  ionized  region  in  the  atmosphere  may  be 
moved  about  by  local  air  currents,  thus  effectively  changing 
the  point  in  space  at  which  the  potential  is  measured.  How- 
ever, it  has  been  well  established  that,  for  measurements 
in  the  outside  atmosphere,  the  fluctuations  can  vary  greatly 
from  time  to  time  (see  Fig.  9).  They  do  not  correlate  with 
gross  wind  velocity.  This  suggests  that  the  composition  of 
the  air  is  involved  and  that  the  nature  of  the  ionic  concen- 
tration in  the  high  conductivity  region  may  affect  the  fluc- 
tuations. 


f Three  experimental  ways  to  reduce  the  fluctuations 

| have  been  attempted:  the  sensor  was  ventilated  by  blowing 

I an  air  stream  across  it,  the  radioactivity  of  the  sensor 
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was  reduced,  the  effect  of  source  symmetry  about  the  hori- 
zontal axis  was  investigated. 

In  general,  laboratory  measurements  where  the 
radioactive  sensor  was  placed  between  capacitor  plates 
have  shown  that  ventilating  the  sensor  with  an  air  stream 
velocity  of  1 to  3 m/s  is  effective  in  reducing  the  fluc- 
tuations. Figure  15  shows  the  effect  of  ventilat.ng  the  sen- 
sor with  an  air  velocity  of  1.5  m/s  (about3mi/h)  origin- 
ating in  a squirrel  cage  blower  1.  85  meters  away..  For  the 
data  in  Fig.  15,  one  sensor  is  grounded  to  the  meter  case 
and  the  other  placed  in  the  center  of  the  capacitor  plates. 
Thus  common  mode  and  dif' .rential  outputs  are  redundant, 
except  for  a factor  of  2.  Unfortunately,  these  results  are 
not  reproducible  day  by  day,  even  in  the  air-conditioned 
laboratory.  Figure  16  shews  data  taken  under  identical 
conditions  4 days  later  on  a very  humid  day.  This  figure 
reveals  that  the  unventilated  sensor  signals  have  greater 
but  slower  fluctuations,  and  that  the  air  velocity  must  be 
increased  to  achieve  a degree  of  smoothing  comparable  to 
that  shown  in  Fig.  15.  One  of  the  significant  results  of 
ventilation  most  easily  observed  in  the  laboratory  is  that 
the  frequency  of  fluctuations  is  always  increased  by  the 
air  flow. 


The  results  of  ventilation  of  the  sensor  in  the  out- 
side atmosphere  are  qualitatively  similar  to  those  in  the 
laboratory.  In  general,  ventilation  reduces  the  amplitude 
and  increases  the  frequency  of  the  fluctuations..  However, 
ventilation  tends  to  be  less  effective  in  the  atmosphere.  On 
rare  occasions  it  do=j  not  help  at  all,  at  least  for  air  ve- 
locities up  to  10  m/s  (approximately  22  mi/h).  Such  an 
occasion  is  shov/n  in  Fig.  17,  where  ventilation  at  10  m/s 
does  not  decrease  the  fluctuation.  These  data,  like  those 
in  the  laboratory,  were  recorded  with  one  sensor  grounded, 
and  the  potential  values  are  distorted  by  the  air  duct  in  the 
proximity  of  the  sensor.  A simultaneous  field  mill  read- 
ing of  the  vertical  field  is  also  shown.  The  data  were  re- 
corded on  a clear,  sunny  day  with  a relative  humidity  of 
43%  and  very  little  wind. 
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The  500  MCi  sensor  that  has  been  used  for  the  above 
measurement  is  much  larger  than  necessary  for  the  high 
impedance  voltmeter.  Measurements  both  in  the  laboratory 
and  in  the  outside  atmosphere  show  that  the  a source 
strength  can  be  reduced  to  about  10  MCi  without  affecting 
the  accuracy  of  the  potential  reading  or  the  time  response 
of  the  instrument.  However,  a reduction  of  a source 
strength  below  this  value  affects  the  accuracy  of  the  po- 
tential reading  and  a source  strength  of  4 MCi  will  produce 
an  error  of  about  5%.  The  reduction  in  source  strength 
does  decrease  the  fluctuation,  but  the  reduction  is  not  dra- 
matic until  a low  enough  value  of  source  strength  is  reached 
to  affect  the  accuracy  of  the  potential  reading.  Figure  18 
shows  the  effect  of  reducing  the  source  strength  to  10  HCi. 
The  decrease  in  fluctuations  is  significant  and  the  error  in 
the  potential  reading  is  negligible.  Under  those  atmospheric 
conditions  (which  generally  prevail)  where  ventilation  re- 
duces the  fluctuations  of  the  500  (*Ci  sensor,  it  also  re- 
duces the  fluctuations  of  the  smaller  sensor. 

As  Ref.  9 points  out,  there  is  an  inherent  error  in 
the  potential  reading  using  the  sensor  shown  in  Fig.  5.  A 
sensor  with  cylindrical  symmetry  about  the  horizontal  axis 
has  been  tried  on  the  high  impedance  voltmeter.  The  dif- 
ference in  the  reading  of  the  atmospheric  potential  is  small 
(anproximately  1 to  2%  for  a sensor  height  of  1 meter), 
and  no  detectable  decrease  in  fluctuations  was  observed. 

It  has  been  found  that  a more  significant  error  in  the 
value  of  the  atmospheric  potential  can  be  caused  by  the 
charging  of  the  fiber  glass  support  rod.  In  the  presence  of 
the  normal  atmospheric  electric  field  (pointed  toward  the 
earth),  this  rod  tends  to  become  negatively  charged . This 
can  produce  an  error  of  several  percent  in  the  potential 
reading.  This  charge  builds  up  to  a stable  value  in  several 
minutes  when  the  electric  field  is  applied  and  requires  sev- 
eral minutes  to  dissipate  after  the  field  is  removed.  If  a 
metal  conductor  2 inches  long  is  placed  around  the  fiber 
glass  arm  and  connected  to  the  guard  circuit  of  the  volt- 
meter (as  shown  in  Fig.  5),  this  charging  effect  disappears. 


.i 

■ i 


I 

1 


u 


60  - 


FIELD  MEASUREMENTS  IN  THE  NEIGHBORHOOD  OF 
OROGRAPHIC  PROTRUSIONS 


The  ground-level  measurements  of  the  atmospheric 
field  in  the  neighborhood  of  orographic  protrusions  were 
performed  using  the  high  impedance  voltmeter  coupled  to 
the  atmosphere  with  radioactive  sensors  and  mounted  on  the 
support  stand  shown  in  Fig.  7.  From  the  results  given  in 
Section  4.  it  is  obvious  that  the  spurious  fluctuations  char- 
acteristic of  this  instrument  will  seriously  degrade  the  ac- 
curacy of  the  measurements,  particularly  of  the  horizontal 
gradient.  Further,  since  such  a series  of  measurements 
requires  a period  of  an  hour  or  more,  it  is  clear  from  the 
results  above  that  significant  changes  in  the  vertical  field 
can  occur  over  this  period  (see  Fig.  14).  Therefore,  dur- 
ing the  course  of  such  a measurement,  it  was  necessary  to 
monitor  the  vertical  field  with  a separate  voltmeter  remote 
from  the  protrusion.  The  voltmeter  readings  were  re- 
corded by  watching  the  voltmeter  with  a telescope  from 
distances  greater  than  30  feet  from  the  instrument. 

Measurements  were  made  of  the  field  perturbations 
in  the  neighborhood  of  walls,  fences,  and  small  isolated 
buildings  located  on  the  grounds  of  the  APL  laboratories 
in  Howard  County,  Maryland.  All  of  these  measurements 
showed  the  expected  result;  that  is,  as  the  protrusion  was 
approached  on  a line  parallel  to  the  ground,  the  atmospheric 
potential  (as  measured  by  the  common  mode  voltage)  de- 
creased and  the  horizontal  gradient  (as  measured  by  the 
differential  voltage)  increased. 

An  example  of  these  results  is  shown  in  Figs.  19 
and  20.  The  high  impedance  voltmeter  with  the  sensors  in 
a horizontal  plane  was  moved  toward  a steel  mesh  fence, 

3. 2 meters  high,  at  altitudes  of  1 and  2 meters.  The  cal- 
culated curves  are  scaled  from  the  results  of  Fig.  2.  Con- 
sidering the  fluctuations  in  the  meter  readings  (about  ±5 
volts)  and  the  time  variation  of  the  atmospheric  potential 
(which  varied  ±20%  during  the  time  of  this  measurement), 
the  agreement  of  experiment  and  theory  is  reasonably 
good. 
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From  the  results  shown  above,  the  field  meter 
should  be  a superior  instrument  for  detecting  the  pertur- 
bations of  the  atmospheric  field  caused  by  orographic  pro- 
trusions. The  field  mill  that  was  available  was  a heavy  in- 
strument, requiring  400  cycle  AC  power,  and  no  field  mea- 
surements with  it  were  attempted.  However,  it  was 
mounted  on  the  laboratory  roof  and  its  performance  as  an 
obstacle  detector  was  compared  to  that  of  the  high  imped- 
ance voltmeter.  A stepladder  covered  with  a fine  con- 
ducting mesh  was  moved  toward  each  instrument.  Figure 
21  shows  the  results  for  the  high  impedance  voltmeter 
and  Fig.  22  for  the  field  mill.  The  signals  are  plotted  as 
a function  of  time  for  distances  of  infinity,  5,  4,  3,  2,  and 
1 meter  between  the  ladder  and  the  instrument.  At  best, 
for  the  voltmeter,  the  change  in  the  horizontal  gradient  can 
be  detected  3 meters  away;  for  the  field  mill,  changes  in 
the  horizontal  gradient  of  2 V/m  are  easily  detected  at  a 
distance  of  5 meters. 


The  orographic  protrusion  measurements  from  the 
Cessna  337  aircraft  were  performed  by  flying  north  and 
south  along  both  the  eastern  and  western  ridges  of  South 
Mountain  in  western  Maryland.  As  was  pointed  out  in  Sec- 
tion 3,  for  a given  flight  it  was  necessary  to  bias  the  air- 
craft frame  so  that  the  signal  channels  to  be  recorded  were 
not  saturated.  All  the  flight  data  reported  here  were  ac- 
quired on  flight  paths  near  the  top  of  the  ridge  and  at  sev- 
eral hundred  feet  east  or  west  of  the  mountain  slope.  There 
was  no  way  to  determine  the  position  of  the  aircraft  except 
for  the  operator's  verbal  description.  These  descriptions 
were  tape  recorded  and  synchronized  with  the  instrumenta- 
tion tape  recorder. 

An  example  of  the  signals  recorded  in  the  aircraft, 
in  this  case  the  lower  wing  sensor  differential  signal 
V2  - V4  flown  over  level  terrain,  is  shown  in  Fig.  23. 

There  is  no  evidence  of  fast  fluctuations  on  the  graph,  and 
the  difference  signal  is  correlated  with  the  gyro  roll  signal. 
Near  the  end  of  this  recording  the  aircraft  is  entering  a 
right  bank.  The  sense  of  the  signal  Vg  - V4  is  positive  for 
a right  bank. 
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Common  Mode  and  Differential  Voltage  versus  Time,  Measurements 
Made  with  High-Impedance  Voltmeter  at  Various  Distances  in  Front 
of  a Ladder  2.4  Meters  High 
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Fig.  22  Verticil  and  Horizontal  Electric  Field  versus  Time,  Meesurements 
Midi  with  Field  Mill  at  Various  Distances  in  Front  of  a Ladder 
2.4  Meters  High 
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Figure  24  shows  the  lower  wing  sensor  differential 
signal  for  passes  north  and  south  along  the  east  side  of 
South  Mountain.  Traveling  south  the  right  wing  is  at  a 
lower  potential  than  the  left  if  the  field  contours  are  as 
those  of  Fig.  3.  and  therefore  the  differential  signal 
V2  - V4  corresponds  in  sense  to  that  of  a left  bank.  The 
voltage  values  given  in  Figs.  23  and  24  correspond  to  the 
difference  signal  as  read  by  the  high  impedance  voltmeter, 
but  not  precisely  to  the  value  of  the  difference  in  the  atmo- 
spheric potential  between  the  positions  of  sensors  V2  and 
V4  since  the  geometric  form  factors  Of  1 (r)  and  0i(r)  of  Eq. 
(22)  have  not  been  determined.  Thus  these  data  do  not  give 
an  accurate  value  of  the  horizontal  gradient.  However,  for 
purposes  of  comparison  with  the  calculations  in  Fig.  3,  it 
suffices  to  measure  the  differential  signal  V2  - V4  for  a 
known  roll  angle  over  level  terrain  and  to  assume  tliat  the 
vertical  field  does  not  change  significantly  during  the  time 
of  a pass  up  and  down  the  mountain.  It  has  been  well  es- 
tablished by  Hill  and  Whyte  (Ref.  10)  that  for  roll  angles  up 
to  30°  the  differential  signal  V2  - V4  is  linear  to  the  ac- 
curacy of  the  measurement.  In  effect,  this  establishes 
for  data  (such  as  that  in  Fig.  24)  the  slope  of  the  equipoten- 
tial  line,  which  can  then  be  converted  to  a horizontal  gradi- 
ent for  a nominal  value  of  the  vertical  gradient  of  100  V/m. 

Averaging  the  horizontal  gradient  for  flights  north 
and  south  along  the  east  and  west  sides  of  the  mountain, 
the  results  are: 


Altitude 

East  side 

West  side 

1300  feet  (396  meters) 
2000  feet  (610  meters) 

14. 1 V/m 
8.1  V/m 

17.5  V/m 
10.  6 V/m 

Ref.  10.  M.  L.  Hill  and  T.  R.  Whyte,  Investigations 
Related  to  the  Use  of  Atmospheric  Electric  Fields  for 
Aircraft  and  RPV  Stabilization.  APL/JHU  TG  1280  (to  be 
published). 


j 


- 68  - 


THK  JOHN*  HOPKINS  U»MV*l»fTY 

APPLIED  PHYSICS  LABORATORY 

tUVl*  I HIM*  MMUM 


Another  flight  was  made  below  and  off  the  mountain  peak  j 

(about  350  meters  altitude  and  several  hundred  feet  from  • 1 

the  slope)  that  yielded  the  following  results  when  averaged 
over  north  and  south  flight  paths  on  both  sides  of  the  moun-  j 

taint  *' 


; Sensor  position  (see  Fig.  8) 

Horizontal  gradient 

i 

; v2  - v4 

19.  5 V/m 

CO 

> 

1 

> 

19.9  V/m 

\ V5  - v6 

16.8  V/m 

/ v2  • v8 

21.3  V/m 

• » 

Given  the  uncertainty  in  the  aircraft  position,  these  ; ! 

results  are  reasonably  consistent  with  those  to  be  expected  *i  ij 

from  the  potential  contours  of  Fig.  3. 


FLUCTUATIONS  IN  THE  ATMOSPHERIC  FIELD 

The  yearly  and  daily  fluctuations  of  the  atmospheric 
electric  field  as  discussed  in  Section  2 and  Fig.  14  show 
that  significant  short  time  fluctuations  can  occur  over 
periods  of  several  minutes.  Chapter  V of  Ref.  3 discusses 
these  short-time  fluctuations  in  some  detail,  and  Refs.  11 
and  12  discuss  recent  investigations  of  the  short-time  fluc- 
tuations in  the  fair  weather  field  made  without  the  use  of 
radioactive  potential  equalizers.  The  results  of  Refs.  11 
and  12  are  in  general  agreement  with  the  limited  amount  of 
data  that  have  been  obtained  at  APL. 

Ref.  11.  K.  Knott  and  G.  Schumann,  "Short  Period 
Variation  of  the  Atmospheric  Electric  Potential  Gradient,  " 
Arch.  Met.  Geoph.  BiokL.  Series  A,  Vol.  21,  1972, 
pp.  319-327. 

Ref.  12.  D.  G.  Yerg  and  K.  R.  Johnson,  "Short 
Period  Fluctuations  of  the  Fair  Weather  Electric  Field,  " 

J.  Geophys.  Res.,  Vol.  79,  No.  15,  1974,  p.  2177. 
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Figure  25  is  a plot  of  the  vertical  and  horizontal 
electric  fields  recorded  from  2:30  PM  on  27  June  until 
9:00  AM  on  28  June.  The  plots  were  recorded  from  the 
roof  of  an  APL  building  with  the  field  mill  mounted  1 meter 
above  the  roof  and  the  high  impedance  voltmeter  mounted 
3 meters  above  the  roof.  The  high  impedance  voltmeter 
was  slower  than  the  one  shown  in  Fig.  4 and  the  data  were 
recorded  on  a low  precision  recorder.  The  results  track 
well  with  those  of  the  field  mill. 

The  vertical  field  has  several  peaks  around  mid- 
night which  were  great  enough  to  overdrive  the  measuring 
systems.  Also  there  are  field  reversals  at  1:00  and  4:30 
AM.  The  horizontal  field  deviated  from  zero  only  during 
these  radical  changes  in  the  vertical  field.  These  data 
were  recorded  on  a cloudy  night,  and  the  measurements 
were  terminated  when  rain  began  to  fall  about  9:00  AM. 

No  measurable  rain  (<0.01  inch)  fell  during  the  night 
There  is  no  variation  in  the  meteorological  data  that 
matches  the  fluctuations  of  the  vertical  electric  field,  al- 
though it  can  be  observed  that  violent  excursions  in  the 
field  do  correlate  with  high  relative  humidity. 

A more  detailed  recording  of  the  time  dependence  of 
the  vertical  and  horizontal  field  changes  for  a field  rever- 
sal is  shown  in  Fig.  26. 

The  high  impedance  voltmeter  station  on  the  flat 
roof  was  operated  for  several  months  prior  to  August  1974, 
and  the  results  can  be  compared  with  meteorological  data 
(wind  velocity,  wind  direction,  precipitation,  humidity, 
temperature,  and  barometric  pressure)  recorded  at  the 
same  location.  The  voltmeter  recorded  both  the  vertical 
and  horizontal  field,  but  a portion  of  the  horizontal  field 
data  are  vitiated  by  a faulty  recorder.  There  is  a further 
question  of  how  accurate  the  electric  field  data  are  during 
and  following  periods  of  precipitation.  However,  typical 
results  can  be  summarized  below. 
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Mswuremant  of  Vortical  and  Horizontal  Atmospheric  Electric 
swing  Field  Reversal,  Measurements  Made  with  Field  Mill  on 
PL.  Time:  Midnight,  27-28  June  1974 


MJCMHOMMWM 

APPLIED  PHYSICS  LABORATORY 

[mii  ihm  > ..O 


For  the  period  of  15  to  26  March,  there  were  three 
periods  of  field  reversals.  Two  of  these  coincide  with 
periods  of  precipitation  and  one  with  a period  of  very  high 
humidity.  One  period  of  precipitation  was  not  accompanied 
by  a field  reversal.  In  general,  the  horizontal  field  is  zero 
except  for  periods  of  violent  changes  in  the  vertical  field. 
During  the  period  of  16  to  22  July,  one  field  reversal  wa« 
observed,  the  only  period  for  which  disturbances  in  the  hori- 
zontal gradient  were  noted.  The  electrical  field  disturbances 
occurred  near  the  time  of  slight  precipitation.  The  results 
of  an  unusual  night  are  shown  in  Fig.  25. 

Data  recorded  and  analyzed  at  A PL  indicate  that 
when  field  reversals  occur  they  persist  for  less  than  10% 
of  the  total  time.  Estimates  of  atmospheric  field  rever- 
sal times  made  by  R.  Weiss  of  APL  from  American  and 
Russian  data  (Refs.  13  and  14)  indicated  that  no  easy  gen- 
eralization could  be  made  from  such  information.  The 
data  from  Pensacola,  Florida  indicated  that  for  a two- 
year  period  field  reversals  occurred  less  than  5%  of  the 
time.  At  Norfolk,  Virginia  the  field  was  reversed  for  15% 
of  the  time  during  1959  and  5.  7%  of  the  time  during  1960. 
Russian  data  showed  variations  from  7.  8 to  20. 5%. 

The  results  of  Refs.  11  and  12  for  the  fluctuations 
of  the  fair  weather  are  in  general  agreement  with  those 
recorded  at  APL.  At  times,  appreciable  fluctuations  are 
observed  over  times  of  several  minutes  for  the  vertical 


Ref.  13.  E.  M.  Trent  and  R.  V.  Anderson,  "Data 
from  a System  of  Atmospheric  Electricity  Stations,  ' 

U.S.  Naval  Research  Laboratory  Report  6162,  20  May  1965. 

Ref.  14.  "Results  of  Ground  Observations  of  Atmo- 
spheric Electricity,"  January-December  1970,  "USSR 
Chief  Administrations  of  Hydro-Meteorological  Service, 
Leningrade,  1972. 
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Held;  there  are  periods  of  hours  when  the  fluctuations  do 
not  exceed  ±10%.  As  is  pointed  out  in  Kef.  11,  fluctuations 
in  the  horizontal  field  are  rare.  The  natural  assumption 
is  that  these  fluctuations  are  of  local  origin  and  are  caused 
by  localized  variations  of  charge  density  or  conductivity  in 
space  floating  over  the  field  sensors  (see  Refs.  11  and  12). 
Since  changes  in  the  horizontal  gradient  are  not  generally 
observed  from  ground-based  measurements,  the  fluctua- 
tions must  be  assigned  to  changes  in  local  conductivity  or 
to  space  charged  clouds  at  some  height  above  the  appara- 
tus. It  must  be  pointed  out  that  the  only  sensitive  evidence 
for  the  high  degree  of  accuracy  with  which  the  horizontal 
gradient  maintains  its  zero  value  over  flat  terrain  is  based 
upon  measurements  made  within  a few  meters  of  the  earth's 
surface.  This  does  not  preclude  the  possibility  that  at  alti- 
tudes of  a hundred  meters  or  so  appreciable  horizontal  gra- 
dients might  be  detected. 
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5.  CONCLUSIONS 


The  results  of  the  Terrain  Avoidance  Program  can 
be  summarized  as  follows: 

1.  The  character  and  magnitude  of  the  atmo- 
spheric electric  fields  about  orographic  pro- 
trusions as  computed  from  Poisson's  equation 
have  been  experimentally  verified  for  a limited 
range  of  parameters. 

2.  Computations  of  the  fields  about  orographic 
protrusions  show  that  both  the  horizontal  and 
vertical  electrical  fields  contain  information 
of  value  for  terrain  avoidance. 

3.  The  solution  of  Poisson's  equation  for  the 
atmospheric  electric  field  in  the  neighborhood 
of  South  Mountain  shows  that  the  horizontal 
electric  field  is  about  5 V/m  for  a ridge  200 
meters  high  as  measured  from  a point  1000 
meters  from  the  ridge  at  an  altitude  of  100 
meters.  Since  5 V/m  is  believed  to  be  a rea- 
sonable horizontal  field  value  to  detect,  and 
the  contour  of  South  Mountain  is  typical  of 
mountains  in  the  eastern  United  States,  then 
field  disturbances  for  this  class  of  mountains 
can  be  detected  at  horizontal  distances  of  five 
times  the  height  for  horizontal  flight  paths 
lying  below  the  peak  of  the  ridge. 

4.  Fluctuations  in  the  electric  field  have  been 
found  to  be  the  limiting  factor  in  detecting 
electric  field  distortions  caused  by  orographic 
protrusions.  These  fluctuations  may  be  real 
or  instrumental.  The  instrumental  errors  are 
not  important  for  field  mills  nor  have  they  been 
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shown  to  be  important  for  high  impedance  volt- 
meters coupled  to  the  atmosphere  with  radio- 
active sensors  for  velocities  of  80  to  160  mi/h. 
Experiments  indicate  thai  it  may  be  possible, 
by  the  proper  choice  of  radioactive  source 
activity  and  sensor  ventilating,  to  greatly  re- 
duce or  eliminate  the  spurious  noise  that  is 
measured  with  radioactive  sensors  in  static 
and  low  velocity  measurements. 

5.  It  seems  reasonable  to  assume  that  the  normal 
fair  weather  atmospheric  field  exists  for  at 
least  90%  of  the  time  over  much  of  the  earth's 
surface.  Perturbations  in  the  horizontal  com- 
ponent of  the  field  are  very  rare  near  the 
ground. 

Measurements  of  the  atmospheric  electric  field  re- 
quire specialized  equipment  and  laborious  data  reduction. 
The  program  reported  here  has  been  limited  in  equipment 
and  manpower,  and  the  following  problems  require  further 
investigation. 

1.  Can  horizontal  field  gradients  of  a few  volts 
per  meter  be  reliably  measured  from  aircraft 
at  distances  of  five  scale  heights  from  typical 
eastern  mountains?  This  problem  includes 
such  factors  as  determining  the  possible  var- 
iations in  atmospheric  conductivity  over  ranges 
of  thousands  of  meters,  eliminating  the  effects 
of  wing  and  tail  motion  of  the  aircraft,  and 
using  the  techniques  (which  are  now  available) 
of  computerized  data  reduction. 

2.  Verify  from  an  aircraft  the  predicted  changes 
in  the  vertical  potential  gradient  near  moun- 
tains. 

3.  Determine  from  aircraft  or  missiles  the  most 
accurate  instrumentation  for  detecting  oro- 
graphic protrusions.  The  accuracy  of  the  field 
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mill  is  well  known  but  further  tests  should  be 
made  with  the  cheaper  high  or  low  impedance 
voltme'er  coupled  to  the  atmosphere  with  radio- 
active sensors. 


4.  Determine  for  a practical  system  the  logic  re- 
quired to  process  the  simultaneous  changes  in 
the  vertical  and  horizontal  fields  as  orographic 
protrusions  are  approached. 

5.  Acquire  more  information  on  the  fluctuations  of 
the  atmospheric  field  as  a function  of  weather, 
altitude,  and  geographic  location.  The  required 
instrumentation  (portable  battery  powered  re- 
corders, sky  cameras,  and  high-altitude  plat- 
forms) is  not  presently  available. 

The  final  result  of  such  a program  would  be  a rec- 
ommended list  of  the  instrumentation  required  for  terrain 
avoidance  for  a vehicle  of  specified  dynamic  performance. 

In  conclusion,  it  should  be  pointed  out  that  the  atmo- 
spheric electric  field  may  be  useful  for  the  following  appli- 
cations: 

1.  RPV  stabilization  (this  is  now  well  established 
over  a wide  range  of  conditions), 

2.  Terrain  avoidance  for  low  altitude  missiles, 

3.  Long  time  drift  correction  in  vertical  gyros, 

4.  Helicopter  stabilization, 

5.  High  tensior  power  line  detection,  and 

6.  Rotation  rate  and  angular  attitude  determin- 
ation for  rotating  missiles. 
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APPENDIX  A* 

ATMOSPHERIC  ELECTROSTATIC  FIELD  NEAR  A 
MOUNTAIN  RIDGE 


INTRODUCTION 


In  the  absence  of  various  disturbing  effects  such  as 
a nearby  thunderstorm,  local  air  pollution,  thermal  in- 
version, and  so  forth,  it  is  well  known  that  the  atmospheric 
electrical  conductivity  may  be  expressed  in  the  form 
a = aoe0y  where  o0  and  0 are  constants  and  y is  the  al- 
titude. Although  there  is  a wide  scatter  in  the  data  avail- 
able for  cr0  and  0,  a representative  set  of  values  is 


a =4.55x10  14  m-  mho 
o 

0 = 2. 22  x lo  4 meter  1 . 


(Al) 


The  electrostatic  potential  is  then  determined  by  the 
equations 


V • J = o , 


(A2) 


i 

r 


i i 


I 1 

1 T 


J = "E  , (A3) 

subject  to  the  applicable  boundary  condition,  where  J is 
the  electric  current  density  and  E the  field  intensity. 
Equations  (A2)  and  (A3)  give  the  desired  differential  equa- 
tion 


v 


2 


<p  + v f n a • 


V(p  = o . 


(A4) 


* Appendix  A was  originally  published  as  APL/JHU  internal 
memo  CJvi  -74  -050  by  Kwang  Yu,  dated  30  April  1974. 
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Although  the  ground  conductivity  is  finite,  it  is  known  that 
the  ground  may  be  represented  by  a perfectly  conducting 
boundary. 


TWO  DIMENSIONAL  PROBLEM 


We  assume  that  the  ground  boundary  exhibits  a 
translational  symmetry  in  one  direction.  Taking  the  x- 
axis  normal  to  this  direction  and  y-axis  vertically  upward, 
the  differential  equation  becomes 


^ <*no> 


ay 

&y 


(A5) 


Before  we  proceed  to  solve  this  problem  with  a suitable 
set  of  boundary  conditions,  we  consider  the  Laplace 
equation: 


ay2 


0 , 


<A6) 


j 

j 


subject  to  the  same  boundary  conditions  and  a twice  differ- 
entiable map  F:0  -•  G = Sipce  VG  = F'V^i,  we  have 

v2g  = F' v2 0 + F"  (Vtf)2  = - [ - 7$)  . VG,  i.  e. , 

F' 

V2G+  [-  Yi  v*]  • 7G=  0 . (A7) 

Comparing  Eq.  (A7)  with  Eq.  (A4),  we  note  that  a mapping 
F,  determined  subject  to  the  condition 

F'  = , K = constant  (A8) 
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will  yield  the  solution  of  the  differential  Eq.  (A5).  In  par- 
ticular, for  the  case  of  plane  symmetry,  # <*  y,  and  hence 


dF 

Tr  = constant  e 
d i/i 


-0y 


-Py 


F - Cj  - C2  e 


In  general,  an  exponential  transformation  of  0, 
-M 


F = Cj  - C2  e 


G(x,  y) 


(A9) 


leads  to 


a = constant  e 

Next  suppose  <P  is  a solution  of  Eq.  (A4) 


(A10) 


7 cp  + 7 in  j • 7^  = 0 


Then  for  any  twice  differentiable  map  F:<p-G  - F(<p)  = G(x,  y), 
we  have 


72G  = F1  72<p  + F"  (7 ip)2 


,F" 


= F'(-7ino-  7<p)  + (^r  7(p)  • (F>  7<p) 


#Ko. 


-7G  • 7 in  (■pT’),  K = constant  . 


Thus  under  the  mapping  F,  0 maps  into  a solution  of  another 
Poisson's  equation  withoreplaced  by  Ko/F'  where  K is  a 
constant. 


i 
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We  now  consider  the  effect  of  boundary  condition. 
Let  us  examine  the  behavior  of  <p  - <p(x,  y),  satisfying 
Eq.  (4),  under  the  conformal  map 

f:z(=x  + iy)  -*  C = f(z)  = £ + iV  . 

Denoting 

0(£(x,  y),  V(x,  y))  - <p(x,y) 
r(£(x,y),  n(x,  y))  = <J(x,  y), 
a simple  calculation  shows  that 

0 = v2<p  + Vino  • ViP  * ((  2 + f 2) 

x y 


^+ttnr>€VWnTV»J 


Thus  the  conformal  mapping  £ + iTJ  = f(x  + iy)  will  produce 
another  solution  to  the  Poisson's  equation.  The  boundary 
condition  satisfied  by  the  new  solution  0 is  determined 
jointly  by  that  for  the  solution  </>  and  the  conformal  mapping. 
We  note  that  at  the  transformed  point,  the  conductivity 
remains  the  same. 

The  results  established  so  far  show  that  the  order 
of  applying  a suitable  conformal  mapping  (to  produce  the 
desired  boundary  effect)  and  a twice-differentiable  map  F 
(to  produce  the  effect  of  the  variable  conductivity)  may  b»= 
interchanged  in  obtaining  an  approximate  solution  to  a 
class  of  two  dimensional  electrostatic  problems. 


BOUNDARY  CONDITION  FOR  A MOUNTAIN  RIDGE 


Suppose  0 is  a solution  of  the  Laplace  equation  with 
the  perfectly  conducting  plane  boundary  condition.  Then 
(in  suitable  units)  0 = y,  and  all  the  equipotential  surfaces 
are  characterized  by  y = constant.  As  is  well  known  from 
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the  electrostatics  in  two  dimensions,  conformal  mapping 
may  be  applied  to  the  coordinates  x and  y to  generate  the 
solutions  to  various  boundary  value  problems.  Thus,  under 
the  mapping 

x + iy  -•  C = u(x,  y)  + iv(x,  y)  . 

The  imaginary  part  v(x,  y)  gives  the  solution  to  the 
electrostatic  potential  problem  with  the  boundary  given  by 

v(x,  y)  = 0 , 

and  the  equipotential  surfaces  are  then  determined  by 

v(x,y)  = constant  . 

A procedure  related  to  the  above  is  the  Schwarz- 
Christoffel  transformation.  * As  an  example,  we  will  con- 
sider the  problem  of  determining  the  electric  potential  <P 
with  a perfectly  conducting  plane  boundary  with  a (grounded) 
perfectly  conducting  wall  of  height  "a"  standing  normal  to 
the  x-axis.  The  transformation  we  seek  is  the  one  charac- 
terized by  the  differential  equation 


dz 

dw 


Vr2  2 

w -a 


(z  = x + iy,  w = u + iv)  . 


Then,  in  the  transformed  plane,  the  problem  reduces  to: 

& + *!*»  o 

*U2  bv2 
= 0 at  u = 0 


♦See,  for  example,  P.  M.  Morse  and  H.  F^rshback, 
Methods  of  Theoretical  Physics.  Vol.  1,  McGraw-Hill, 
New  York,  1953,  p.  445. 
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|^-  = constant,  say,  1 at  v = 0 


4 1 = 0 at  v = 0 . 

The  differential  equation  leads  to  w = V t?  + a?.  As  the 
solution  in  w-plane  is  equal  to  the  imaginary  part  v,  the 
potential  in  the  z plane  is  given  by 


4>U,  y)  = Im  < V(*2  + z - y2)  + i 2xy) 


..  \ W1  ! J72  2 2.2  , , 2 2 , 2 2.  2.)  ..... 

0(x,  y)  = Tg  I -y+z)+4xy  - (x  -y+a)|  (All) 

The  equipotential  surfaces  are  given  by  0 = constant. 


Explicitly,  setting  4>  ~ Sc,  we  obtain 


r 

x.  c) = y : 


'4c2  + c(x2  + a2) 

X2  + C 


Note  that  any  of  the  equipotential  surface  may  be  taken  as 
the  ground  boundary  provided  that  <p  is  corrected  by  Sc  to 
satisfy  zero  potential  condition  at  the  ground. 


APPLICATION  TO  SOUTH  MOUNTAIN 

Figure  A1  is  a typical  cross  section  of  South  Moun- 
tain north  of  Harper's  Ferry.  This  profile  can  be  approxi- 
mated very  closely  by  the  choice  c = 0. 56  and  a correspon- 
dence 


x - X = 1270x  (feet) 


w 


,2  Comparison  of  Cross  Section  of  South  Mountain  with  Analytic 
Approximation,  See  Fig.  A1 
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y-  Y = 1270y  - 225  (feet) 

in  Eq.  (A12).  The  analytic  approximation  and  the  actual 
cross  section  as  given  in  Fig.  A1  are  compared  in  Fig. 

A2. 

To  obtain  the  solution  corresponding  to  the  South 
Mountain  region,  we  perform  the  exponential  transfor- 
mation (cf.  Eq.  (A9))  to  Eq.  (All).  The  constants  are  so 
chosen  that  the  electric  field  intensity  at  the  surface  at 
infinity  is  100  V/m.  The  equipotential  surfaces  computed 
are  shown  in  Fig.  A3. 

In  Fig.  A4,  the  surfaces  on  which  the  magnitudes 
of  electric  field  intensity  are  constant  (|2  | constant)  are 
shown  for  the  South  Mountain  region.  In  contrast  to  the 
equipotential  surfaces  of  Fig.  3 which  are  as  expected,  the 
surfaces  of  constant  |2  | at  low  altitudes  exhibit  an  inter- 
esting behavior.  In  particular,  the  surface  for  |E  | = 91.  1 
V/m  contains  the  saddle  point  of  the  function  |E(x,  y)|. 
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Horizontal  Range  (meters) 

Fig.  A3  Equipotential  Surfaces,  E0  * 100  V/m 
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APPENDIX  B* 

ELECTROSTATIC  POTENTIALS  NEAR  A PERIODIC 
SERIES  OF  MOUNTAIN  RIDGES 


In  a previous  memorandum  (Appendix  A)  an  analyti- 
cal technique  is  considered  which  may  be  used  to  construct 
solutions  to  a large  class  of  two  dimensional  atmospheric 
electrostatic  problems  each  with  an  infinitely  conducting 
(not  necessarily  planar)  ground  boundary.  The  method 
consists  of  constructing  a suitable  conformal  transforma- 
tion (including  the  Schwarz -Christoffel)  of  a uniform  half- 
space solution  to  satisfy  the  ground  boundary  condition  and 
then  performing  an  exponential  mapping  to  obtain  the  expo- 
nential variation  of  the  atmospheric  conductivity. 

In  this  memorandum  an  application  of  this  method 
is  presented  for  a class  of  boundary  conditions  corres- 
ponding to  a periodic  series  oi  mountain  ridges  (such  as 
the  Appalachian  or  seme  parts  of  the  Rocky  mountains). 

We  then  consider  the  effects  of  a localized  anomalous  air 
conductivity.  Such  regions  occur  as  a result  of  some  local 
radioactivities  and/or  other  localized  pollutants. 


PERIODIC  SERIES  OF  MOUNTAIN  RIDGES 
Boundary  Condition 

The  first  task  is  the  selection  of  a suitable  moun- 
tain profile  that  is  sufficiently  simple  for  computational 
ease.  The  profile  selected  to  do  this  can  be  obtained  from 
the  imaginary  part  of  the  conformal  mapping 

7.  - fiz)  = z - ai  eli>Z,  o > 0 . (Bl) 


•Appendix  B was  originally  published  as  APL/JRU  internal 
memo  QM-  74-103  by  Kwasag  Yu,  dated  24  July  1974. 
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i.  e. , 


n - y - a cos  bx  e = 0 . (B2) 

Evidently,  in  order  that  this  represents  the  profile  of  a 
periodic  series  of  mountain  chain  the  solution  y * y(x) 
must  exist  and  is  continuous  for  all  x.  This  requirement 
is  clearly  satisfied  if  ab  < e"1  *. 


Exponential  Mapping  and  the  Field  Distribution 

In  order  that  the  exponential  mapping  tj  -»  F(e^) 
generates  the  solution  to  the  atmospheric  electrostatic  po- 
tential, we  put 


an  E„  -foix.y) 

<P  = F(e"^)  = - ~~  (1-e  ),  P > 0 (B3) 

E = constant  . 

o 

<P  * 0 at  the  ground  boundary  is  represented  by  ij  = 0.  The 
atmospheric  conductivity  is  o = ffoexp(+0u);  the  small  per- 
turbation introduced  into  the  atmospheric  conductivity  as 

e-a  cos  bx  eby  damps  0ut  rapidly  as  the  altitude  (y)  in- 
creases. The  constant  /3  then  establishes  the  relation  be- 
tween the  dimensionless  coordinates  (x,y)  and  the  physical 
coordinates  (horizontal  range,  altitude). 


♦This  may  be  seen  from  the  function  X(u,  v)  = ueu  - <*  cos 
v = 0,  a * real.  The  first  term  has  a minimum  at  u - -1, 
and  for  all  u > -1,  it  is  a continuous  monotonically  in- 
creasing function  of  u.  The  second  term  continuously 
oscillates  with  the  amplitude  a . Thus  the  solution  u = u(v) 

exists  only  if  |a|  < e . As  |a|  - e , - ”. 

v = » 


J 

.*1 


. s 

I 1 


1 
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The  electric  field  is  given  by 


E = - V<p  = E ' 7T?(x,y)  , 
o 


and  the  family  of  equipotential  surfaces, 

© = constant  = - F , 
o 

is  obtained  from  the  equation 


’l  = Jfn(E-T¥T)-  ^ 

^ o o 


In  Fig.  Bl,  a number  of  equipotential  surfaces  correspond 
to  several  constant  values  of  V- 


LOCALIZED  ANOMALOUS  AIR  CONDUCTIVITY 

As  the  starting  point,  the  electrostatic  field  is  con- 
sidered over  a perfectly  conducting  boundary  with  the  con- 
stant conductivity  replaced  by 
two  uniform  conductivity  re-  Region  || 

gions.  The  semicircular  in- 
ner region  has  a radius  R.  In  S' 

these  regions,  the  potential  / 1 /q  \ 

distribution  $ is  given  by  the  ///'////  Tr/zTl  7 / 
imaginary  part  of  the  function  t'  ' ' ' ' ' ' ' ' ' ' ' ' ' ' 


|z|  < R (B5) 


f2(z>=  -Czt  |z|>R 
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rv 


zg2 

gl  + ff2  y 


2 2 
x + y s R 


« = 


(B6) 


x2  + y2  > R 

1 2 x + y 


Several  equipotential  surfaces  for  this  potential  distribution 
function  are  shown  in  Fig.  B2  for  the  case  ^ 1 /o' 2 = 2. 

Under  the  conformal  mapping  z -•  z - ai  eibz,  the 
conductivity  boundary  maps  into 

_2  2 . 2 2 -2by  . „ . , . . -by 

R=x+y+ae  J + 2a(x  sin  by-y  cos  dx)  e J = constant. 

and  the  potential  $ maps  into  V 

2o. 


1 • 


g +2g  (y-a  COS  bx  e'by) 


*2  " V2  = '<y‘a  008  bx  e"by>  { vT5"2 


gl-a2 


(B7) 


(x  + a sin  bx  e"bx)  + (y  - a cos  bx  e’by)2  ^ 


We  now  consider  the  exponential  map 


<P  = F(e'^>  = - — (1  - e'?V)  . 


(B8) 


Note  that  the  equipotential  surfaces  of  <p  and  V coincide  as 
before.  In  Fig.  B3  several  equipotential  surfaces  are 
shown  which  correspond  to  those  given  in  Fig.  B2.  In  the 
region  where  the  curves  are  dotted  the  potentia*  varies 
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Fig.  B3  Equipotantial  Surfaces  over  a Periodic  Series  of  Mountain  Ridges  with  a 
Region  of  Anomalous  Conductivity,  o^/cj  = 2,  a = 1,  b = ir/10,  R = 2 
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rapidly.  This  is  due  t<-  the  discontinuous  conductivity  that 
was  used  and  will  be  a:  ..  -t  if  the  conductivity  is  made  to 
vary  continuously. 


SURFACES  OF  CONSTANT  |E| 

As  was  done  in  Ref.  Bl,  the  geometry  of  the  sur- 
faces on  which  |e|  remains  constant  is  examined.  In 
Fig.  B4  three  such  surfaces  are  shown  for  /3  = 1 (for  ease 
of  computation).  This  amounts  to  setting  the  distance 
scale  so  that  unit  distance  corresponds  to  450C  meters. 

To  examine  the  shapes  of  the  "equi-|E|"  surfaces 
with  respect  to  those  of  equipotential  surfaces,  the  sine* 
of  the  angle  0 between  |E|  = constant  surface;  the  equipo- 
tential surface  is  calculated  along  the  equipotential  sur- 
faces corresponding  to  V - 0,  1,  2,  3,  and  4 with  0 set 
equal  to  1.  The  results  are  shown  in  Figs.  B5  through  B7. 
To  examine  the  0 dependence,  for  the  same  set  of  values 
for  »},  sin  6 was  computed  with  0 = 1/8.  The  results  are 
shown  in  Figs.  B8  through  BIO.  In  Figs.  B5  to  BIO,  the 
solid  lines  show  the  results  corresponding  to  the  presence 
of  anomalous  conductivity  (o  ^ /<r 2 = 2),  and  the  dotted  lir.es 
to  the  normal  conductivity  <crj /cr2  = 1). 

As  expected,  for  increasing  V the  amplitude  of 
|sin  ©|  decreases,  while  the  maximum  amplitude  increases 
as  0 decreases  from  1 to  1/8.  The  rapid  variation  near 


♦The  sine  of  the  angle  © between  the  two  surfaces  is  com 
puted  from 


e • (ExVE2) 

Sin  © » — x-  , 

|E|  • | VE| 

where  E*  = E • E and  ez  is  the  right-handed  unit  vector 
normal  to  the  xy  plane. 


Ref.  Bl.  J.  A.  Stratton,  Electromagnetic  Theory. 
McGraw-Hill,  1941. 
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the  conductivity  boundary  will  probably  disappear  if  the 
conductivity  is  made  to  vary  smoothly. 

To  investigate  the  behavior  of  |e|  as  a function  of 
altitude  in  the  valley  region  we  put  Enx  = 10)  = X. 

X = E2(x  = 10,  y.  0) 

= e-20(y+ae-by)(1_abe-by)2  . 

At  the  ground,  W(x  - 10)  = 0 = y + ae"by,  and  hence 

= *2<1  + byJ  W1  + byJ2  + b2yJ  . 

* V * 0 

where  y0  denotes  the  ground  value  of  y at  x = 0.  Setting 

b2y0 

V (r+~byVw°-792,  (a=  b = ff/10> 


we  have 


/ > 0 if  0 < | 

I « 0 if  0 * | 


^x=10  y = y 

* 10>  y y0  (<  o if  0 > 0c . 

In  case  0 is  less  than  0C,  the  amplitude  of  |e|  is  an  in- 
creasing function  of  altitude  at  the  ground  level  in  the 
valley.  To  see  the  behavior  off  the  ground  we  examine 


(|i) 

8y  x = 10,  y > y_ 


-20(y  + ae  by) 


(1  - ab  e"by)  • [0(1  - ab  e~by) 


/ 


m 
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Since  ab  < e , the  above  vanishes  only  if  the  square 
bracket  term  vanishes.  This  occurs  at  y = yc  given  by: 


bv  = in 


2/3  + b-Kb  + lb/3 


For  y > y . grl  < 0 . 

c y vx  = 10 

Again  assuming  /3  > /3C,  we  consider  the  behavior 
of  |E|  at  x = 10,  y = yc.  Evaluating  X^and  Xyy,  we 
have 


2 -byc  -2/Jnc 

2ab  e e c 


-ab  e n + b 


x = 10.  y = y 


x = 10,  y = y 


x = 10,  y = y„ 


> 0 , 


-byc  -23tJc  2 -byc  r -by  -i 

-2(l-abe  )e  ab  e 2/3(l-ab  e c)  + b|  < 0 


where  V = t/l 

'c  "x  = 10,  y = y 


showing  that  the  point  (x  =10,  y = yc)  is  a saddle  point. 


Values  of  Parameters  for  Appalachian  Mountains 

For  the  Appalachian  mountains  in  West  Virginia,  a 
typical  peak-to-peak  distance  is  on  the  order  of  5 to  10 
miles  (8  to  16  km'  with  a typical  peak-to-valley  height  in 
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the  range  of  1000  to  3000  feet.  Taking  the  height  to  "wave- 
length" ratio  to  be  1/10  we  see  that  the  surface  correspond- 
ing to  V - 0.  5,  a = 1,  b = ff  / 10  may  be  used  to  represent 
the  ground  boundary.  The  value  of  P appropriate  for  this 
case  is  # = 0.  111. 


Further  Application  of  the  Analytic  Techniq 


Recent  measurements  of  the  atmospheric  electro- 
static field  both  at  APL  (near  ground)  and  in  the  vicinity  of 
South  Mountain  (airborne  measurements)  indicate  a fluc- 
tuation in  E of  as  much  as  10  V/m.  This  may  be  due  to 
locally  fluctuating  regions  of  anomalous  conductivity 
caused,  for  example,  by  some  local  concentrations  of  air 
pollutants. 

A two-dimensional  analytical  analog  of  this  may  be 
constructed  by  considering  a circular  (in  two-dimensions) 

region  of  anomalous  conduc-  

tivity.  The  solution  corres-  / Anomalous 

ponding  to  this  can  then  be  ( ) p°n^‘nctivity 

transformed  exponentially,  V J 6910,1 

and  the  field  plots  (Ex,  Ey,  ^ ' Normal 

|E|,  etc.)  from  this  can  be 
used  to  correlate  with  the 

experimental  data.  /// 7777777777  Gmund 
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APPENDIX  C* 


MEASUREMENT  OF  A rM^Sr^ERIC  POTENTIAL  IN  THE 
NEIGHBORHOOD  OF  CONDUCTORS  OF  J1MFLE  SHAFE 
WHICH  BEAR  A FREE  CHARGE 


INTRODUCTION 


The  electrostatic  autopilot  proposed  for  ' P'/'e 
measures  the  atmospheric  potential  or  a quantity  proportion- 
al to  it,  relative  to  the  body  of  the  vehicle,  at  the  position 
of  the  sensor.  Since  aircraft  are  known  to  become  charged 
to  potentials  much  higher  or  lower  than  their  surroundings, 
it  is  necessary  to  examine  what  effect  these  potentials  have 
on  the  attitude  sensors. 


It  would  seem  clear  (at  least  for  shapes  that  obey 
certain  symmetry  properties)  that,  if  the  sensors  are 
placed  in  positions  where  they  are  on  equipotential  surfaces 
that  result  from  a charge  on  the  vehicle,  the  difference 
signals  from  the  sensors  will  be  independent  of  the  poten- 
tial of  the  vehicle  relative  to  its  surroundings  and  will  ac- 
curately measure  the  atmospheric  potential.  The  following 
calculations  for  the  sphere  and  prolate  spheroid  show  that 
this  is  true.  However,  these  calculations  also  indicate 
that,  in  practice,  it  may  be  difficult  to  establish  these 
equipotential  surfaces  for  a two-axis  system  and  that  the 
charge  on  the  vehicle  can  either  paralyze  or  destroy  the 
electronics,  even  though  arranged  to  only  measure  the 
atmospheric  potential. 


♦Appendix  C was  originally  published  as  APL/JHU  internal 
memo  -73-167  by  C.  S.  Leffel,  Jr.,  dated  2 November 
1973. 
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CONDUCTING  SPHERE  IN  A UNIFORM  ELECTRIC  FIELD 

For  a conducting  sphere  of  radius  r0  in  a uniform 
electric  field  (see  Fig.  Cl) 

ip  = E x = E r cos  0 , (Cl) 

0 0 o 

the  potential  field  is  given  by 


E r r 

<p  « Er  cos  0 - — = — cos  0 + <p  — 
£ s r* 

r 


where  the  first  term  is  the  given  potential  field,  the  second 
term  is  the  perturbation  term  caused  by  the  charges  in- 
duced on  the  conductor,  and  the  third  term  is  the  potential 
caused  by  any  free  charge  on  the  conductor. 


s 4 w ( r 


The  electric  field  for  r * r0  is  giver,  by 


_ r r 

E=-V<p=(E  cos  0 - 2E  ~-+(P  -?)? 

o o 3 s 2 

r r 


+ (E  sin  0 - E — y sin  6)  0 . 
O O 3 

r 


Along  the  x axis,  the  field  "augmentation  factor,  " neglect- 
ing  ip s,  is  three. 

If  potential  probes  are  placed  along  any  line  through 
the  center  of  the  sphere,  at  equal  distance  r > r0  from  the 
center,  and  if  Vj  = <p( r,  0)  - <pB  and  V2  = <P(r,  6 + v)-<PB, 
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the  voltages  measured  with  voltmeters  inside  the  sphere, 
then 


Vl  + V2  rn 

V = ’s -«• 


and  the  common  mode  voltage  is  a measure  of  the  potential 
of  the  sphere  relative  to  its  surroundings.  The  difference 
of  Vj  and  V2  is 


V!  ‘ V2  = 2Eor  cos  Sii-  “%■>  = 2<po(1*  1 (C6) 


Thus  (Vj  - V2)  measures  ip^  apd  the  measurement  is 
more  sensitive  the  greater  r exceeds  rD.  On  the  other 
hand.  Vj  for  example,  is  given  by 


Vi  = ^0(1--T)+  Vf*1*  * 


and  if  <PB  is  great  enough  and  r is  great  enough,  Vj  ■»  -<Pe, 
which  will  paralyze  or  destroy  any  electronic  system. 

If  the  sensors  are  not  symmetrically  placed,  a 
change  in  <PS  can  produce  a false  error  signal.  For  two 
sensors  along  the  y1  axis,  fixed  to  the  sphere,  placed  at 
distances  rj  and  r2  from  the  center,  then  the  error  signal 
to  maintain  the  y1  axis  horizontal  is,  with  a = (tr/2 > - 0. 


V2  = Eosin  alrj  + rg  - r*  + ~T > * Vo  - “)  • <C8) 


As  long  as  9S  - 0,  the  error  signal  will  drive  the  angle  a 
to  zero,  but  if  the  vehicle  becomes  cliarged,  then  a signal 
Vj  * ^2  be  produced  which  will  drive  the  y1  axis  from 
its  horizontal  position. 
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The  value  of  the  potential  along  the  x axis  for  a 
sphere  of  1 meter  radius  in  a field  of  100  V/m  and  with  a 
surface  potential  of +100  volts  is  as  shown  in  Fig.  Cl  and 
plotted  in  Fig.  C2. 


CONDUCTING  PROLATE  SPHEROID  IN  A UNIFORM 
ELECTRIC  FIELD 

Laplace’s  equation  can  be  integrated  in  terms  of 
elementary  functions  for  the  prolate  spheroid  oriented  in 
an  electric  field  as  shown  in  Fig.  C3.  ( and  £ define  a 
family  of  confocal  ellipsoidal  and  hyperboloidal  surfaces 
rotated  about  the  x axis.  Here  the  potential  of  the  field  of 
constant  E0  is  given  by 

E 

<p  =Ex  = ±-2  <?+  a2)  <C  + b2)  . (C9) 

o o c 


Assume  that  the  spheroid  of  constant  potential  is  defined 
by  f = 0.  The  potential  field  caused  by  free  charge  on  the 
ellipsoid  is  given  by 


(£)  = z~~ 

s 8ffc 


/ 


JL 


(i+  b ) 


2'  V777 


which  when  integrated  yields 
V <£) 


8ire  c in 


1 . Vt  + a2  + c 
+ a2  - c 


(CIO) 


The  potential  field  caused  by  the  charge  induced  on  the 
spheroid  by  the  5 field  is  given  by 
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Fig.  C2  Potential  FMd  along  tha  x-Axis  for  a Sphare,  & = 100  Vota.^j  = -100  V/m  x 


1 


I 


TWC  JOHN*  MOM  MM  UMtVCMlTV 

APPLIED  PHYSICS  LABORATORY 

•*.**»  «MWi  HMTUM 


s?+l;  b?+f 


-b2f>-a2 


x = ±^y«  + a2)(r+a2) 


y = ± 1 >/-({  + b2)  (f  + b2) 

c 


Fig.  C3  Syrtam  of  Coordinate  for  Prolata  Spheroid 
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At  4 = 0,  V>i  (0  ‘ -<PQ  and  for  the  region  exterior  to  the 
spheroid,  the  potential  is  given  by 


= <*>„  + v6<e)  + «s.<4)  . 


The  electric  field  in  the  4 direction  is  given  by 
F?  - .1  - -2  yle  + a2  + b2  dip 

VTT 


hj  34 


34 


and  if  we  set 


k- -l  + IinaJLo 

a c a - c 


the  value  of  the  4 component  of  the  E field  is 


jLi! 

K Vc  + *2 

1 EoV«+a2 

/T7 

( 2 

K 2 

-2 


VTT 


Sb&ib 


«-J~W fcaA  -H|  ^~ff1 
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tt+a2)-c2Jj  <C13) 

Evaluating  Eq.  (C13)  at  x = a,  y = 0 which  means  4=0, 

C = -b2,  the  result  is 

_2c2e 

E(x  = a,  y = 0)  = — j 2°  a+c  * ' <CU) 

-2b  c + ab  in  2 — - 
a - c 

If  a s>  b,  the  field  augmentation  factor  can  be  very  large  at 
the  sharp  end  of  the  ellipsoid. 

The  results  of  Eq.  (C12)  are  laborious  to  plot,  but 
except  for  the  lack  of  spherical  symmetry  are  similar  to 
those  of  the  sphere.  The  first  result  to  note  is  that  the 
equipotential  surfaces,  as  given  by  Eq.  (CIO)  for  the  free 
charge  on  the  conducting  ellipsoid  defined  by  4 = 0,  are 
simply  the  ellipsoids  defined  by  ?>  0.  Since  C = -b2  at 
the  poles  (x,  0)  and  C = -a2  at  the  equator  (0,  y),  from  the 
equations  for  x and  y given  in  Fig.  C3, 


: = * V«+  ^ 


at  the  poles 


y = ± V*+b  at  the  equator. 

Thus,  from  Eq.  (CIO),  if  a = 5,  b = 1,  and  (q/8ire)  = 1.  the 
following  results  are  obtained: 


L 

PB<0 

X 

Z 

0 

0.9358 

5 

1 

l 

0.7984 

5.099 

1.414 

2 

0.7196 

5. 196 

1.032 

4 

0.6229 

5.  385 

2.  236 

6 

0.5614 

5.  567 

2.645 

i 
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The  equipotential  surfaces,  which  must  approach  a sphere 
as  | increases,  lie  further  from  the  conducting  surface 
along  the  y-axis  than  along  the  x-axis. 


The  values  of  the  potential  near  the  ellipsoid  along 
the  x axis  are  plotted  on  Fig.  C4  and  along  the  y axis  on 
Fig.  C5,  for  a * 5 meters,  b = 1 meter,  E0  = 100  V/m 
and  <ps(£  = 0)  = 93.  5 volts.  Compared  to  the  sphere,  the 
rate  of  change  of  potential  is  much  faster  along  the  x axis 
and  slower  along  the  y axis.  In  either  case,  at  distances 
far  enough  from  the  origin,  the  value  of  ip  asymptotically 


approaches  <po. 


If  sensors  are  mounted  along  a line  through  the 
center  of  the  ellipsoid  and  if  V is  measured  relative  to 
<ps(£  = 0),  and  if  Vi  = V(x,y),  V2  = V (-x,  -y),  then  since 
<Po<*.  y)  = - <Po  (-*.  -y).  then  from  Eq.  (C12) 


V1  = <po+  *8®  + = 0) 


V_  = - <p  + <f>(0  - <P.(£)  - Oft  = 0),  and  again  we 

td  OS  L S 


have  the  results 


Vl+V2 


= <Pstt)  - <Pgft  = 0) 


(C15) 

(C16) 


Vj  - v2  = 2[«poce,  0 + c>]  . 

From  Eq.  Cll,  it  is  observed  that  Ojft,  0 can  be  writtenas 


<P.ft.  0 = - <po« . Or© 


Thus  Eq.  (Cl  6)  can  be  written  as 

V1  - v2  = 2<Po«,0  {1  -yft)]  , 


(C17) 
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where  y(|)  is  1 for  £ = 0.  Equations  (C15)  and  (C'17)  are 
identical  in  form  to  Eqs.  (C5)  and  (C6)  for  the  sphere. 
Thus  ips (4),  caused  by  the  self  charge  and  <p0,  caused  by 
the  external  field  can  be  measured  if  the  geometrical  fac- 
tors are  known. 


CONCLUSIONS 

1.  Since  the  potential  field  caused  by  free 
charge  on  the  surface  of  the  vehicle  is  fixed  in  the  coordi- 
nate systems  of  the  vehicle  itself,  it  can  always  be  mea- 
sured by  an  equation  of  the  form  of  Eq.  (Cl 5),  no  matter 
how  complicated,  provided  the  sensors  are  on  equipoten- 
tial  surfaces.  Further,  for  the  difference  voltages,  the 
effect  of  the  free  charge  potential  will  always  cancel  out. 

However,  in  the  case  of  an  actual  RPV,  the 
equipotential  surfaces  caused  ty  free  charge  may  be  very 
difficult  to  determine.  For  the  roll  axis,  sensors  on  the 
wing  tips  will  by  the  symmetry  of  the  vehicle  body  easily 
be  placed  on  or  nearly  on  an  equipotential  surface  cauFed 
by  free  charge.  But  the  problem  is  much  more  difficult 
for  the  pitch  axis:  if  the  fore  and  aft  sensors  are  not  on 
equipotential  surfaces,  an  accumulation  of  charge  on  the 
vehicle  body  will  produce  a climb  or  dive.  To  find  the 
free  charge  equipotential  surfaces  by  analysis  for  a real- 
istic vehicle  snape  seems  a formidable  task.  An  experi- 
mental approach  in  which  the  vehicle  is  suspended  and  the 
sensor  positions  adjusted  for  zero  difference  signals  as 
voltages  are  applied  to  the  vehicle  skin  may  be  more 
promising.  However,  the  problem  is  complicated  by  the 
ground  plane  effects  (image  charges)  which  have  not  been 
considered  here,  and  which  may  be  important  when  an  RPV 
is  charged  under  laboratory  conditions. 

Since  it  is  necessary,  in  order  to  detect  the 
atmospheric  potential,  that  the  sensors  be  placed  outside 
the  free  charge  equipotential  surface  of  the  vehicle,  it  is 
always  possible  that  the  vehicle  can  become  charged  to  a 
potential  that  will  saturate  any  realistic  measuring  circuit. 


> 
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2.  The  above  calculations  for  the  prolate  spher- 
oid are  not  adequate  for  the  case  in  which  the  major  axis  of 
the  spheroid  is  displaced  from  the  direction  of  the  atmo- 
spheric field.  That  is  to  say,  the  values  calculated  for  <pQ 
and  y will  no  longer  be  correct.  However,  because  of  the 
symmetry  of  the  spheroid  the  form  Eq.  (C17)  will  still  be 
applicable.  However,  ewen  in  the  case  of  lack  of  symmetry, 
it  would  appear  that  the  error  signal  Vi  - V2  will  always  be 
such  to  move  the  body  so  that  the  sensors  lie  in  atmospher- 
ic equipotential  planes.  For  in  general,  the  difference 
voltage  will  be  of  the  form: 


V - V = 
v2  V1 


pol(1  -’'l)+<po2(1 


and  since  y < 1 and  = <fio2  = 0 when  the  sensors  are  in 
atmospheric  equipotential  planes,  the  error  signal  can  be 
coupled  to  the  controls  80  that  V2  - Vj  is  driven  to  zero. 
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